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Abstract 

The purpose of this review is to provide basic ingredients of holographic QCD to non- 
, experts in string theory and to summarize its interesting achievements in nuclear and hadron 

,_^ physics. We focus on results from a less stringy bottom-up approach and review a stringy 

f^ top-down model with some calculational details. 
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1 Introduction 

The approaches based on the Anti de Sitter/confornial field theory (AdS/CFT) correspon- 
dence [l}l3] find many interesting possibihties to explore strongly interacting systems. The 
discovery of D-branes in string theory [4] was a crucial ingredient to put the correspondence on 
a firm footing. Typical examples of the strongly interacting systems are dense baryonic matter, 
stable/unstable nuclei, strongly interacting quark gluon plasma, and condensed matter systems. 
The morale is to introduce an additional space, which roughly corresponds to the energy scale 
of 4D boundary field theory, and try to construct a 5D holographic dual model that captures 
certain non-perturbative aspects of strongly coupled field theory, which are highly non-trivial 
to analyze in conventional quantum field theory based on perturbative techniques. There are in 
general two different routes to modeling holographic dual of quantum chromodynamics (QCD). 
One way is a top-down approach based on stringy D-brane configurations. The other way is 
so-called a bottom-up approach to a holographic, in which a 5D holographic dual is constructed 
from QCD. Despite the fact that this bottom-up approach is somewhat ad hoc, it reflects some 
important features of the gauge/gravity duality and is rather successful in describing properties 
of hadrons. However, we should keep in mind that a usual simple, tree-level analysis in the holo- 
graphic dual model, both top-down and bottom-up, is capturing the leading Nc contributions, 
and we are bound to suffer from sub- leading corrections. 

The goal of this review is twofold. First, we will assemble results mostly from simple bottom- 
up models in nuclear and hadron physics. Surely we cannot have them all here. We will devote 
to selected physical quantities discussed in the bottom-up model. The selection of the topics is 
based on authors' personal bias. Second, we present some basic materials that might be useful 
to understand some aspects of AdS/CFT and D-brane models. We will focus on the role of the 
AdS/CFT in low energy QCD. Although the correspondence between QCD and gravity theory 
is not known, we can obtain much insights on QCD by the gauge/gravity duality. 

We organize this review as follows. Section ^ reviews the gauge/gravity. Section ^ briefly 
discuss developments of holographic QCD and demonstrates how to build up a bottom-up model 
using the AdS/CFT dictionary. After discussing the gauge/gravity duality and modeling in the 
bottom-up approach, we proceed with selected physical quantities. In each section, we show 
results mostly from the bottom-up approach and list some from the top-down model. Section [4] 
deals with vacuum condensates of QCD in holographic QCD. We will mainly discuss the gluon 
condensate and the quark-gluon mixed condensate. Section [5] collects some results on hadron 
spectroscopy and form factors from the bottom-up model. Contents are glueballs, light mesons, 
heavy quarkonium, and hadron form- factors. Section [6] is about QCD at flnite temperature and 
density. We consider QCD phase transition and dense matter. Section [7] is devoted to some 
general remarks on holographic QCD and to list a few topics that are not discussed properly in 
this article. Due to our limited knowledge, we are not able to cover all interesting works done in 
holographic QCD. To compensate this defect partially, we will list some recent review articles 
on holographic QCD. 

In Appendix, we look back on some basic materials that might be useful for non-experts in 
string theory to work in holographic QCD. |Aj we review the relation between the bulk mass 
and boundary operator dimension, p] we present a D3/D7 model and axial U(l) symmetry in 
the model. O we discuss non-Abelian chiral symmetry based on D4/D8/D8 model. PmEI we 
describe how to calculate the Hawking temperature of an AdS black hole. |FJ we encapsulate 
the Hawking-Page transition and sketch how to calculate Polyakov loop expectation value in 
thermal AdS and AdS black hole. 

We close this section with a cautionary remark. Though it is tempting to argue that holo- 
graphic QCD is dual to real QCD, what we mean by QCD here might be mostly QCD-like or a 
cousin of QCD. 



2 Introduction to the AdS/CFT correspondence 

The AdS/CFT correspondence, first suggested by Maldacena [l], is a duality between gravity 
theory in anti de Sitter space (AdS) background and conformal field theory (CFT). The original 
conjecture states that there is a correspondence between a weakly coupled gravity theory (type 
IIB string theory) on AdS^ x S^ and the strongly coupled A/" = 4 supersymmetric Yang-Mills 
theory on the four- dimensional boundary of AdS^. The strings reside in a higher-dimensional 
curved spacetime and there exists some well-defined mapping between the objects in the gravity 
side and and the dual objects in the four-dimensional gauge theory. Thus, the conjecture allows 
the use of non-perturbative methods for strongly coupled theory through its gravity dual. 

2.1 Dp brane dynamics 

The duality emerges from a careful consideration of the D-brane dynamics. A Dp brane sweeps 
out {p + 1) world-volume in spacetime. Introducing D branes gives open string modes whose 
endpoints lie on the D branes and the open string spectrum consists of a finite number of 
massless modes and also an infinite tower of massive modes. The open string end points can 
move only in the parallel (p + 1) directions of the brane, see Figure 1(a) , and a Dp brane can 



be seen as a point along its transverse directions. The dynamics of the T)p brane is described 
by the Dirac-Born-Infeld (DBI) action |5j and Chern-Simons term. 



Sd, = -Tp j dP+'x e- V-det (PMafc + 27ra'F,b) + Scs 



(2.1) 



with a dilaton e '^. Here gab is the induced metric on Dp. P denotes the pullback and Fab is 
the world-volume field strength. Tp is the tension of the brane which has the form 
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and it is the mass per unit spatial volume. Here gs is the string coupling and Ig is the string 
length, a' is the Regge slope parameter and related to the string length scale as Ig = vc/. In 
general, states in the closed string spectrum contain a finite number of massless modes and an 
infinite tower of massive modes with masses of order rus = Ij^ = a'~^' ^. Thus, at low energies 
E <C rus, the higher order corrections come in powers of a'E'^ from integrating out the massive 
string modes. If there are a Nc stack of multiple D branes, the open strings between different 



branes give a non-Abelian U{Nc) gauge group, see Figure 1(b) In the low energy limit, we can 



integrate out the massive modes to obtain non-Abelian gauge theory of the massless fields. 
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(a) D3 branes sweep (3 + 1) dimensions 
in (9+1) space time. 




(b) Nc — S stack of D3 branes and all the possible classes 
of open strings. 



Figure 1: The configurations of Nc stack of D3 branes in Wd spacetime. 



Now, we take p = 3 and consider Nc D3 brane stacks in type IIB theory. The low energy 
effective action of this configuration gives a non-Abelian gauge theory with U{Nc) gauge group. 



In addition, this gauge group can be factorized into U{Nc) = SU{Nc) x U{1) and the U{1) part, 
which describes the center of mass motion of the D3 branes, can be decoupled by the global 
translational invariance. The remaining subgroup SU{Nc) describes the dynamics of branes 
from each other. Therefore we see that in the low energy limit, the massless open string modes 
on Nc stacks of D3-branes constitute AA = 4 SU{Nc) Yang-Mills theory I6| with 16 supercharges 
in (3 + 1) spacetime. From (2.1) for p = 3, we obtain the effective Lagrangian at low energies 
up to the two derivative order 



C 



1 



'i-^gs 



-Tr 



-FauF^"" + -D 



^l^i 



ijjt^^i 



1 -•• wi2 , ^T,/ 



[0^</'^■ 



+ -^'T^'D^^j 



-¥r[(l)i,^i]] (2.3) 



with a gauge field Af^, six scalar fields </>* and four Weyl fermions ^^. 

In fact, the original system also contains closed string states. The higher order derivative 



corrections for the Lagrangian (2.3) come both in powers of a'E from the massive modes and 



powers of the string coupling gsNc for loop corrections. It is known that the string coupling 
constant gs is related by the 10-dimensional gravity constant as G^^^' ~ g'^l^ and thus the 
dimensionless string coupling is of order G^^^'E^, which is negligible in the low energy limit. 
Therefore, at low energies closed strings are decoupled from open strings and the physics on the 
Nc D3 branes is described by the massless A/" = 4 super Yang-Mills theory with gauge group 
SU{Nc). 

2.2 AdS^ X S^ geometry 

Now we view the same system from a different angle. Since D branes are massive and carry 
energy and Ramond-Ramond (RR) charge, Nc D3 branes deforms the spacetime around them to 
make a curved geometry. Note that the total mass of a D3 brane is infinite because it occupies 
the infinite world-volume of its transverse directions, but the tension, or the mass per unit 
three-volume of the D3 brane 

n = l/{2nfgslt (2.4) 

is finite. 

In the flat spacetime, the circumference of the circle surrounding an origin at a distance r 
is 27rr, and it simply shrinks to zero if one approaches the origin. But if there is a stack of 
D3 branes, it deforms the spacetime and makes throat geometry along its transverse directions. 
Thus, near the D3 branes, the radius of a circle around the stack approaches a constant R, an 



asymptotical infinite cylinder structure, or AdS^ x S^, see Figure 2(b) The Nc D3 brane stack 
is located at the infinite end of the throat and this infinite end is called the "horizon" . In the 
near horizon geometry, a D3 brane is surrounded by a five-dimensional sphere S^. 
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(a) Flat spacetime for r ":$> R. (b) Nc D3 branes deform the spacetime. 

Figure 2: The two descriptions of the Nc D3 configuration. 

To be more specific, let us start with type IIB string theory for p = 3. We find a black hole 
type solution which is carrying charges with respect to the RR four-form potential. The theory 



has magnetically charged D3 branes, which are electrically charged under the potential dA^ and 
it is self-dual *Fr, = F^. The low energy effective action is 

S = (^ / d'^xV^g ^e-'HR + <V^f) " |i^l) • (2-5) 

We assume that the metric is spherically symmetric in seven-dimensions with the RR source at 
the origin, then the Nc parameter appears in terms of the five-form field RR-field strength on 
the five-sphere as 

/ *F5 = N, (2.6) 

where S^ is the five-sphere surrounding the source for a four- form field C4. Now by using the 
Euclidean symmetry we get the curved metric solution |7||9l for the D3 brane 



where 



m = 1 + § (2i 



with the radius of the horizon R 



R^ = ^ATTQsNca' = ^AttqsNcII . (2.9) 

dn^ is the five-sphere metric. For r ^ R we have /(r) ~ 1 and the spacetime becomes flat 
with a small correction R'^/r'^ = AngsNclg/r'^. This factor can be interpreted as a gravitational 
potential since G^^^) ~ g^l^ and Md3 ~ N^T^ ~ Njoslt thus R'^/r^ ~ GMos/r'^- In the near 
horizon limit, this gravitational effect become strong and the metric changes into 

ds^ = ^r^f'^dx^dx^ + ^ {dr^ + r^d^l) . (2.10) 

This is AdS5 X S^. 

The geometry by the D3 branes is sketched in Figure |2(b)[ Far away from the D3 brane 
stacks, the spacetime is flat (9-|- 1) dimensional Minkowski spacetime and the only modes which 
survive in the low energy limit are the massless closed string (graviton) multiplets, and they 
decouple from each other due to weak interactions. On the other hand, close to the D3 branes, 
the geometry takes the form AdS^ x S^ and the whole tower of massive modes exists there. This 
is because the excitations seen from an observer at infinity are close to the horizon and a closed 
string mode in a throat should go over a gravitational potential to meet the asymptotic fiat 
region. Therefore, as we focus the lower energy limit, the excitation modes should be originated 
deeper in the throat, and then they decouple from the ones in the flat region. Thus in the low 
energy limit the interacting sector lives in AdS^ x S^ geometry. 

2.3 The gauge/gravity duality 

So far, we have considered two seemingly different descriptions of the Nc D3 brane configuration. 
As we mentioned, each of the D3 branes carries the gravitational degrees of freedom in terms of 



its tension, or the string coupling gs as in (2.4). So the strength of the gravity effect due to Nc 



stacks of D3 branes depends on the parameter ggNc- 



If gsNc <C 1, from (2.9) we see that R <^ Is and therefore the throat geometry effect is 
less than string length scale. Thus the spacetime is nearly fiat and the fluctuations of the D3 
branes are described by open string states. In this regime the string coupling gs is small and 
the closed strings are decoupled from the open strings. Here the closed string description is 
inapplicable since one needs to know about the geometry below the string length scale. If we 



take the low energy limit, the effective theory, which describes the open string modes, is A/" = 4 
super Yang-Mills theory with SU{Nc) gauge group. 

On the other hand, if ggNc ^ 1, then the back-reaction of the branes on the background 
becomes important and spacetime will be curved. In this limit the closed string description 
reduces to classical gravity which is supergravity theory in the near horizon geometry. Here the 
open string description is not feasible because gsNc is related with the loop corrections and one 
has to deal with the strongly coupled open strings. Again, if we take the low energy limit, the 
interaction is described by the type IIB string theory in the near-horizon geometry, AdS^ x S^. 

The gauge/gravity correspondence is nothing but the conjecture connecting these two de- 
scriptions of Nc D3 branes in the low energy limit. It is a duality between the AA = 4 super 
Yang-Mills theory with gauge group SU{Nc) and the type IIB closed string theory in AdS^ x S^, 
see Figure [Sl The relation between the Yang-Mills coupling gyM ^^^ ^^^ string coupling strength 
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Figure 3: The sketch of the AdS/CFT correspondence. 



gs is given by 
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Then, the 't Hooft coupling A = gy^^c can be expressed in terms of the string length scale 



A = R^l 



A /,4 
s • 



(2.12) 



Therefore, the dependence of gs^c becomes the question of whether the 't Hooft coupling is 
large or small, or the gauge theory is strongly or weakly coupled. 

The two descriptions can be viewed as two extremes of r. For the sake of convenience, we 



use the coordinate z = E? /r. Then, the AdS^ x S^ metric (2.10) becomes 

R^ 



ds' 



-y {'q^"'dx^dx^ + dz^) + R^d^l, 



(2.13) 



which shows the conformal equivalence between AdS^ and flat spacetime more clearly. In (2.13), 
each z-slice of AdS^ is isometric to four- dimensional Minkowski spacetime. In this coordinate, 
z = IS the boundary of AdS^, where Yang-Mills theory lives, with identifying x^ as the 
coordinates of the gauge theory. If z — >• oo, the determinant of the metric goes to zero and it 
is the Poincare horizon. Here the factor E? / z"^ also has some relation with the energy scales. 
If the gauge theory side has a certain energy scale E, the corresponding energy in the gravity 
side is {z/R)E. In order words, a gauge theory object with an energy scale E is involved with 
a bulk side one localized in the z-direction at 2: ~ 1/E [T, 10, 11 . Therefore, the UV or high 
energy limit corresponds to z — )• (or r — )• 00) and the IR or low energy limit corresponds to 
z — )• 00 (or r — )• 0). 

The operator-field correspondence between operators in the four-dimensional gauge theory 
and corresponding dual fields in the gravity side was given in |[2ii3j. Then the AdS/CFT corre- 
spondence can be stated as follows, 

/j.Jd^xM'-^)Oix)\ =Zsugra (2.14) 

where 4>oix) = (j){x,u — t- oo) and the string theory partition function Zgugra at the boundary 
specified by (pQ has the form 

(2.15) 



•^sugra 



The relation (2.14) implies that the generating functional of gauge-invariant operators in CFT 



can be matched with the generating functional for tree diagrams in supergravity. 

3 Holographic QCD 

Ever since the advent of the AdS/CFT correspondence, there have been many efforts, based 
on the correspondence, to study non-perturbative physics of strongly coupled gauge theories in 
general and QCD in particular. 

Witten proposed [12] that we can extend the correspondence to non-supersymmetric theories 
by considering the AdS black hole and showed that this supergravity treatment qualitatively well 
describes strong coupled QCD (or QCD-like) at finite temperature: for instance, the area law 
behavior of Wilson loops, confinement/deconfinement transition of pure gauge theory through 
the Hawking-Page transition, and the mass gap for glueball states. In [13] symmetry breaking by 
expectation values of scalar fields were analyzed in the context of the AdS/CFT correspondence, 
which is essential to encode the spontaneous breaking of chiral symmetry in a holographic QCD 
model. Regular supergravity backgrounds with less supersymmetries corresponding to dual 



confining M = 1 super Yang-Mills theories were proposed in (14 15 . It has been shown by 
Polchinski and Strassler [16] that the scaling of high energy QCD scattering amplitudes can be 
obtained from a gravity dual description in a sliced AdS geometry whose IR cutoff is determined 
by the mass of the lightest glueball. Important progress towards flavor physics of QCD has been 
made by adding flavor degrees of freedom in the fundamental representation of a gauge group 



to the gravity dual description 17 . Chiral symmetry breaking and meson spectra were studied 
in a nonsupersymmetric gravity model dual to large Nc nonsupersymmetric gauge theories [18], 
where flavor quarks are introduced by a D7-brane probe on deformed AdS backgrounds. Using 



a D4/D6 brane conflguration, the authors of 19 explored the meson phenomenology of large N, 



QCD together with U(1)a chiral symmetry breaking. They showed that the chiral condensate 
scales as l/niq for large niq. A remarkable observation made in 19 is that in addition to the 



confinement/deconfinement phase transition the model exhibits a possibility that another tran- 
sition set by Tfund could happen in deconflned phase, T > Tdeconf, where Tdeconf = Mkr/C^t^)- 
Since the value of Mkk is around 1 GeV, we can estimate Tfi^^onf ~ 160 MeV. In this case for 
Tdeconf < T < Tfund there exist free unbound quarks and meson bound sates of heavy quarks 
and above Tf^nd the meson states dissociate into free quarks, which in some sense mimics the 



dissociation of heavy quarkonium in quark-gluon plasma (QGP). However, we should note that 
meson bound states in Dp/Dq systems are deeply bound, while the heavy quarkonia in QCD 
are shallow bound states. In this sense the bound state that disappears above Tfy^^d could be 
that of strange quarks rather than charmonium or bottomonium 1 20| . 

To attain a realistic gravity dual description of (large Nc) QCD, non-Abelian chiral sym- 
metry is an essential ingredient together with confinement. Holographic QCD models, which 
are equipped with the correct structure for the problem, namely, chiral symmetry and con- 
finement, have suggested in top-down and bottom-up approaches. They found to be rather 
successful for various hadronic observables and for certain processes dominated by large Nc- 
Based on a D4/D8/D8 model, Sakai and Sugimoto studied hadron phenomenology in the chiral 



limit ruq = 0, and the chiral symmetry breaking geometrically 121 22 . More phenomenolog- 



ical holographic QCD models were proposed 23-25 . In |23[|24|, chiral symmetry breaking is 



realized by a non-zero chiral condensate whose value is fitted to meson data from experiments. 
Hadronic spectra and light-front wave functions were studied in |26| based on the "Light-Front 
Holography" which maps amplitudes in extra dimension to a Lorentz invariant impact separa- 
tion variable C in Minkowski space at fixed light-front time. Light-Front Holography has led to 
many successful applications in hadron physics including light-quark hadron spectra, meson and 
baryon form factors, the nonperturbative QCD coupling, light-front wave-functions, see p7] - |29| 
for a review on this topic. In [30], a relation between a bottom-up holographic QCD model and 
QCD sum rules was analyzed. 

Now, we demonstrate how to construct a bottom-up holographic QCD model by looking at 
a low-energy QCD. For illustration purposes, we compare our approach with the (gauged) linear 
sigma model. The D3/D7 model is summarized in Appendix [B| with some calculational details. 



For a review of the linear sigma model, we refer to 31 . Some material in this section is taken 
from |32|. Suppose we are interested in two-flavor QCD at low energy, roughly below 1 GeV. 
In this regime usually we resort to the effective models or theories of QCD for analytic studies 
since the QCD lagrangian does not help much. 

To construct the holographic QCD model dual to two flavor low-energy QCD with chiral 
symmetry, we first choose relevant fields. To do this, we consider composites of quark fields that 
have the same quantum numbers with the hadrons of interest. For instance, in the linear sigma 
model we introduce pion-like and sigma-like fields: vf ~ qr'j^q and a ~ qq, where r is the Pauli 
matrix for isospin. In the AdS/CFT dictionary, this procedure may be dubbed operator/field 
correspondence: one-to-one mapping between gauge-invariant local operators in gauge theory 
and bulk fields in gravity sides. Then we introduce 

qL-f''t''qL^Al^{x,z), 

rRqi^{2/z)X-P{x,z). (3.1) 

An interesting point here is that the 5D mass of the bulk field is not a free parameter of the 
model. This bulk mass is determined by the dimension A and spin p of the dual 4D operator in 
AdSrf+i. For instance, consider a bulk field X(x, z) dual to q{x)q{x). The bulk mass of X{x^ z) 
is given by m?^ = ( A — p) ( A + p — d) with A = 3, p = and d = 4, and so m^ = — 3. For more 
details, see Appendix \A\ 

To write down the Lagrangian of the linear sigma model, we consider (global) chiral symmetry 
of QCD. Since the mass of light quark ~ 10 MeV is negligible compared to the QCD scale 
^QCD ~ 200 MeV, we may consider the exact chiral symmetry of QCD and treat quark mass 

effect in a perturbative way. Under the axial transformation, q — )• e"*'^^^ ' ^g, the pion-like and 
sigma-like states transform as vf — )■ vf + 0cj and a ^ a — 6 ■ n. From this, we can obtain terms 
that respect chiral symmetry such as vf^ + cr'^. Similarly we ask the holographic QCD model 
to respect chiral symmetry of QCD. In AdS/CFT, however, a global symmetry in gauge theory 



corresponds local symmetry in the bulk, and therefore the corresponding holographic QCD 
model should posses local chiral symmetry. This way vector and axial- vector fields naturally fit 
into chiral Lagrangian in the bulk as the gauge boson of the local chiral symmetry. 

We keep the chiral symmetry in the Lagrangian since it will be spontaneously broken. Then 
we should ask how to realize the spontaneous chiral symmetry breaking. In the linear sigma 



model, we have a potential term like ((vf +0" 



^2\2 



that leads to spontaneous chiral symmetry 



breaking due to a nonzero vacuum expectation value of the scalar field o^ (a) = c. In this case 
the explicit chiral symmetry due to the small quark mass could be mimicked by adding a term 
—ea to the potential which induces a finite mass of the pion, m^ ~ e/c. In a holographic QCD 
model, the chiral symmetry breaking is encoded in the vacuum expectation value of a bulk scalar 

rUgZ + (z^, 



23 



24 



it is given by (X) 



field dual to qq. For instance in the hard wall model 
where rUg and (^ are proportional to the quark mass and the chiral condensate in QCD. In the 
D3/D7 model, chiral symmetry breaking can be realized by the embedding solution as shown in 
Appendix |B] 

The last step to get to the gravity dual to two flavor low-energy QCD is to ensure the 
confinement to have discrete spectra for hadrons. The simplest way to realize it might be to 
truncate the extra dimension at z = Zm such that the radial direction z of dual gravity runs 
from zero to Zm- Since the radial direction corresponds to an energy scale of a boundary gauge 
theory, 1/zm maps to Aqcd- 

Putting things together, we could arrive at the following bulk Lagrangian with local SU(2)lx 



SU(2)r, the hard wall model 23,24 



^HW 



d X dz^ Tr 



'45i 



(F| + F|) + |DX|2 + 3|X| 



(3.2) 



where Dy^X = d^,X - iAl^ X + IXAr^ and Al,r = Al^^t" with Tr(t"t'') = ^5"^. The bulk 
scalar field is defined by X = Xq g^*'^ * , where Xq = {X). Here 55 is the five dimensional gauge 
coupling, (75 = 12-k'^/Nc. The background is given by 



ds^ 



1 



{dr - dx^ - dz^) , 0<z<z„ 



(3.3) 



Instead of the sharp IR cutoff in the hard wall mode, we may introduce a bulk potential that 
plays a role of a smooth cutoff. In [33], this smooth cutoff is introduced by a factor e~ with 
^{z) = z^ in the bulk action, the soft wall model. The form ^i^z) = z^ in the AdS would ensure 
the Regge-like behavior of the mass spectrum m^ ~ n. The action is given by 



5'sw 



d X dze ^Jg Tr 



'4^1 



(F| + F|) + |DX|2 + 3|X| 



(3.4) 



Here we briefly show how to obtain the 4D vector meson mass in the soft wall model. The 
vector fleld is deflned hy V = Al + Ar. With the Kaluza-Klein decomposition V^{x,z) = 
95 Y.n Vn{z)p1^{x), we obtain 



B„ 



dzie dzVn) + rn^e Vn = 0, 



(3.5) 



where B = ^{z) — A{z) = z +log2; in the AdS geometry (3.3). With Vn = e ipn, we transform 



the equation of motion into the form of a Schrodinger equation 



(3.6) 



where V{z) = z^ + 3/(4z^). Here m„ is the mass of the vector resonances, and p meson 
corresponds to n = 0. The solution is well-known in quantum mechanics and the eigenvalue m^ 



is given by 33 



m„ 



4c(n + l), 



(3.7) 



where c is introduced to restore the energy dimension. 

The finite temperature could be neatly introduced by a black hole in AdS^+i, where d is the 
dimension of the boundary gauge theory. The background is given by 



1 / „. X ..9 ,^2 dz 



^* = ^ Ui^)dt' - dx' - — 1 , (3.8) 

where f{z) = 1 — z*^ / zf^. The temperature of the boundary gauge theory is identified with the 
Hawking temperature of the black hole T = d/{A'KZh)- In Appendix O and ^ we try to explain 
in a comprehensive manner how to calculate the Hawking temperature of a black hole. 

Now we move on to dense matter. According to the AdS/CFT dictionary, a chemical poten- 
tial in boundary gauge theory is encoded in the boundary value of the time component of the 
bulk U(l) gauge field. To be more specific on this, we first consider the chemical potential term 
in gauge theory, 

^t^ = f^qQU- (3.9) 

Then, we introduce a bulk U(l) gauge field A^ which is dual to g7^(?. According to the dictionary, 
Aq{z — )• 0) ~ ciz'^~'^~P + C2Z^~P, we have Ao{z — )• 0) ~ iJ,q. In the hard wall model, the solution 
of the bulk U(l) vector field is given by 

At{z) = fi + pz^ , (3.10) 

where p and p are related to quark chemical potential and quark (or baryon) number density in 
boundary gauge theory. It is interesting to notice that in chiral perturbation theory, a chemical 
potential is introduced as the time component of a gauge field by promoting the global chiral 



symmetry to a local gauge one 34 



4 Vacuum structures 

At low energy or momentum scales roughly smaller than 1 GeV, r > 1 fm, QCD exhibits 
confinement and a non-trivial vacuum structure with condensates of quarks and gluons. In this 
section, we discuss the gluon condensate and quark-gluon mixed condensate. 

The gluon condensate {Ct^^Ga^) was first introduced, at zero temperature, in J35| as a mea- 
sure for nonperturbative physics in QCD. The gluon condensate characterizes the scale symmetry 
breaking of massless QCD at quantum level. Under the infinitesimal scale transformation 

^). = (1 + '^^)^M ' 
A'^ = {l-5X)A^, 

ci = {l-^-5\)q, (4.1) 

the trace of the energy momentum tensor reads schematically 

d,j^ = r; ^ -{^G%Gn . (4.2) 

Here J^ is the dilatation current, as is the gauge coupling, T^y is the energy- momentum tensor 
of QCD. Due to Lorentz invariance, we can write (T^^) = —e^acV^^ ■, where e^ac is the energy 
of the QCD vacuum. Therefore, the value of the gluon condensate sets the scale of the QCD 
vacuum energy. In addition, the gluon condensate is important in the QCD sum rule analysis 
since it enters in the operator product expansion (OPE) of the hadronic correlators [35 1. At high 
temperature, the gluon condensate is useful to study the nonperturbative nature of the QGP. 
For instance, lattice QCD results on the gluon condensate at finite temperature l36j indicate that 
the value of the gluon condensate shows a drastic change around Tc regardless of the number 
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of quark flavors. The change in the gluon condensate could lead to a dropping of the heavy 
quarkonium mass around Tc [37| . 

In holographic QCD, the gluon condensate figures in a dilaton profile according to the 
AdS/CFT since the dilaton is dual to the scalar gluon operator Tt{G^uG^^). The 5D grav- 
ity action with the dilaton is given by 



5 = 7 



1 

2^ 



d^x^ 



n+^-'^dM^d'^'cP 



(4.3) 



where 7 = +1 for Minkowski metric, and 7 = —1 for Euclidean signature. We work with 
Minkowski metric for most cases in this paper. The solution of this system is discovered in [38l|39| 
by solving the coupled dilaton equation of motion and the Einstein equation: 



ds^ 






(4.4) 



and the corresponding dilaton profile is given by 



(^(z) 



log 



1 + cz^ 



1 



cz^ 



+ 



(4.5) 



where (/)q is a constant. At 2; = 1/c^'^ there exists a naked singularity that might be resolved in 
a full string theory consideration. Near the boundary 2; — )• 0, 



{z) ~ cz^ 



(4.6) 



Therefore, c is nothing but the gluon condensate up to a constant. Unfortunately, however, c 
is an integration constant of the coupled dilaton equation of motion and the Einstein equation 
and therefore, it will be determined by matching with physical observables. In |39| , the value 
of the gluon condensate is estimated by the glueball mass. An interesting idea based on the 
circular Wilson loop calculation in gravity side is proposed to calculate the value of the gluon 
condensate G2 = (^G^^G';^^) [iol. The value is determined to be G2 = 0.010 ± 0.0023 GeV 



at zero temperature 40 . A phenomenological estimation of the gluon condensate in QCD sum 



rules gives (^G;j^Gr) ^ 0.012 GeV 



35 



Now we consider quark-gluon mixedcondensate {qaf^^G'^'^q), which can be regarded as an 
additional order parameter for the spontaneous chiral symmetry breaking since the quark chi- 
rality flips via the quark-gluon operator. Thus, it is naturally expressed in terms of the quark 
condensate as 



{qa^uG'^'^q) = ml{qq) ■ 



(4.7) 



In 



41 



an extended hard wall model is proposed to calculate the value of rnl. The bulk action 



of the extended model is given by 



S = I d^Xy/gTi 



|L>XP + 3|XP 



4ff; 



■(F| + F^) + |Z?$|2-5$2 



(4.8) 



where <^ is a bulk scalar field dual to the 4D operator on the left-hand-side of Eq. (4.7). Then 



the chiral condensate and the mixed condensate are encoded in the vacuum expectation value 
of the two scalar fields. 



{X{x,z)) = -{mz + az^), 
{^{x,z)) = -{ciz-^ +aMZ^) 

D 



(4.9) 
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where ci is the source term for the mixed condensate and gm represents the mixed condensate 
(^M = {<iR(^ fivG^^ Ql) ■ Taking ci = 0, a source-free condition to study only spontaneous 
symmetry breaking, we determine the value of the mixed condensate or rriQ by considering 
various hadronic observables. In this sense the mixed condensate is not calculated but fitted 
to experimental data like the chiral condensate in the hard wall model. The favored value of 
the rriQ in [4l| is 0.72 GeV^. A new method to estimate the value of m^ in is suggested in 



42 



where a nonperturbative gauge invariant correlator (the nonlocal condensate) is calculated in 
dual gravity description to obtain rrig. With inputs from the slopes for the Regge trajectory 
of vector mesons and the linear term of the Cornell potential, they obtained ttiq = 0.70 GeV^, 
which is comparable to that from the QCD sum rules, 0.8 GeV^ 



43 



5 Spectroscopy and form factors 

Any newly proposed models or theories in physics are bound to confront experimental data, 
for instance hadron masses, decay constants and form factors. In this section, we consider the 
spectroscopy of the glueball, light meson, heavy quarkonium, and hadron form factors in hard 
wall model, soft wall model, and their variants. 

5.1 Glueballs 

Glueballs are made up of gluons with no constituent quarks in them. The glueball states are in 
general mixed with conventional qq states, so in experiments we may observe these mixed states 



only. Their existence was expected from the early days of QCD 44,45 . For theoretical and 
experimental status of glueballs, we refer to (46|[47| . 

The spectrum of glueballs is one of the earliest QCD quantities calculated based on the 
AdS/CFT duality. In [12] , Witten confirmed that the existence of the mass gap in the dilaton 
equation of motion on a black hole background, implying a discrete glueball spectrum with 
a finite gap. Extensive studies on the glueball spectrum were done in [48j|49j and there also 
comparisons between the supergravity results and lattice gauge theory results were made. 

Now we consider a scalar glueball (O"*"*") on R x S^ as an example 48 . When the radius 



of the circle S^ is very small i? — t- 0, only the gauge degrees of freedom remains and the gauge 
theory is effectively the same as pure QCD3 Jil2||48j. Using the operator/field correspondence, 
we first find operators that have the quantum numbers with glueball states of interest and then 
introduce a corresponding bulk field to obtain the glueball masses. In this case we are to solve 
an equation of motion for a bulk scalar field 0, which is dual to tri^^ in the AdSs Euclidean 
black hole background. The equation of motion for (j) is given by 



5^(V55.</>0 = 0, (5.1) 



and the metric is 



b^Y\2 (2 b^ 

7) ''^V ? 



ds^=(p^-^) dp^+(p^-^)dT^ + p^ + {dxf + d^l , (5.2) 



where r is for the conipactified imaginary time direction. For simplicity, we assume that <J3 is 



independent of r 12,48 and seek a solution of the form (j){p,x) = f{p)e , where k is the 



momentum in R^. Then the equation of motion for /(p) reads 



P dp^P ^^^%'^'^ ^^' ^^-^^ 



where m? is the three-dimensional glueball mass, m^ = —k'^ 12,48]. By solving this eigenvalue 



equation with suitable boundary conditions, regularity at the horizon (p = b) and normalizability 
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f ^ p ^ at the boundary {p — )• oo), we can obtain discrete eigenvalues, the three-dimensional 
glueball masses. In the context of a sliced AdS background of the Polchinski and Strassler set 



up 16 , which is dual to confining gauge theory, the mass ratios of glueballs are studied in 50 



More realistic or phenomenology-oriented approaches follow the earlier developments. In the 
soft wall model the mass spectra of scalar and vector glueballs, and their dependence on the 
bulk geometry and the shape of the soft wall are studied in [51] . The exact glueball correlators 
are calculated in (52] , where the decay constants as well as the mass spectrum of the glueball 
are also obtained in both hard wall and soft wall models. Here we briefly summarize the scalar 
glueball properties in the soft wall model |51 , 52 . Following a standard path to construct a 



bottom-up mode, we introduce a massless bulk scalar field (j) dual to the scalar gluon operator 



Ti{F^yF^^^) to write down the bulk action as 51 



S ^ 



Jd'^e- 



'9MN{dMmd''<P) 



(5.4) 



where <!• = z^ as in the soft wall model. The equation of motion for (p(q, z) can be transformed 
to a one dimensional Schrodinger form 



< - V{z)tlj 



q^ip, 



(5.5) 



where ip = e' 



-{<S>+3lnz)/2. 



with g^ 



m 



51 



The glueball mass spectrum is then given as the 
eigenvalue of the Schrodinger type equation with regular eigenfunction at z = and z = oo 

2 



m„ 



4(n + 2)c, 



(5.6) 



where n is an integer. 



n 



0,1,2, 



c is introduced to make the exponent $ dimensionless. 



$ 



cz"^. 



and it will be fit to hadronic data. Since the vector meson mass in the soft wall model 
4(n + l)c, we calculate the ratio of the lightest (n=0) scalar glueball mass w-gg to the 

The properties of the glueball at finite temperature 



53 



51 



and also in the soft wall model |53 54 by calculating 



IS m^^ 

p meson mass to obtain mQQ/m'i = 

is studied in the hard wall model 

the spectral function of the glueball in the AdS black hole backgroudn. The spectral function 

is related to various Green functions, and it can be defined by the two-point retarded Green 

function as p{uj,q) = —2lTa.G{uj,q). The retarded function can be computed in the real-time 

AdS/CFT, following the prescription proposed in 55 . Both studies using the soft wall model 



predicted that the dissociation temperature of scalar glueballs is far below the deconfinement 



(Hawking-Page) transition temperature of the soft wall model. See Section 6.1 and in Appendix 



[F| for more on the Hawking-Page transition. Note that below the Hawking-Page transition 
temperature, the AdS black hole is unstable. In [53], the melting temperature of the scalar 
glueball from the spectral functions is about 40 — 60 MeV, while the deconfinement temperature 
of the soft wall model is about 190 MeV t56j. This implies that we have to build a more refined 
holographic QCD model to have a realistic melting temperature 53 , 54 . 



5.2 Light mesons 

There have been an armful of works in holographic QCD that studied light meson spectroscopy. 
Here we will try to summarize results from the hard wall model, soft wall model and their 
variants. 

In Table [T| we list some hadronic observables from hard wall models to see if the results 
are stable against some deformation of the model. In the table, where * means input data and 
the model with no * is a fit to all seven observables: Model A and Model B from the hard wall 



model 1 23] , Model I from a hard wall model in a deformed AdS geometry 57 , and Model II from 
a hard wall model with the quark-gluon mixed condensate 41 .In 24 , the following deformed 
AdS background is considered 



ds' 



vr 



2ZmSm['Kz/{2Zm)] 



{dt^ — dxidx^ — dz^) , < z < Zm, 



(5.7) 
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and it is stated that the correction from the deformation is less than 10%. The backreaction 
on the AdS metric due to quark mass and chiral condensate is investigated in l57j . One of the 
deformed backgrounds obtained in |57] phenomenologically reads 



ds"^ = ^e-^^^"^) (dt^ - dxidx' - dz"^) ,0 < z < Zm , (5.8) 



Z2 



where B{z) = -^z'^ + ^^^z^ + fj^^- In Table hi we quote some results from this deformed 
background. Dynamical (back-reacted) holographic QCD model with area-law confinement and 
linear Regge trajectories was developed in [58j. 





Model I 


Model II 


Model A 


Model B 


Experiment 


rrip 


775.8* 


775.8 


775.8* 


832 


775.49 ±0.34 


rriai 


1348 


1244 


1363 


1220 


1230 ± 40 


u 


92.4* 


80.5 


92.4* 


84.0 


92.4 ±0.35 


Fr 


334 


330 


329 


353 


345 ±8 


Fti' 


481 


459 


486 


440 


433 ± 13 


m^ 


139.6* 


139.3 


139.6* 


141 


139.57 ±0.00035 


QpTTTT 


4.46 


4.87 


4.48 


5.29 


6.03 ±0.07 



Table 1: Meson spectroscopy from the hard-wall model and from its variations: Model I 57 , 
Model II [41] , Model A I23l, Model B l23l. The experimental data listed in the last column are 
taken from the particle data group [59|. All results are given in units of MeV except for the 
condensate and the ratio of two condensates. 

We remark that the sensitivity of calculated hadronic observables to the details of the hard 
wall model was studied in [60] by varying the infrared boundary conditions, the 5D gauge 
coupling, scaling dimension of qq operator. It turns out that predicted hadronic observables are 
not sensitive to varying scaling dimension of qq operator, while they are rather sensitive to the 
IR boundary conditions and the 5D gauge coupling [60| . 

In addition to mesons, baryons were also studied in the hard wall model (6l| - [64| . It is pointed 
out in 62 , 63 that one has to use the same IR cutoff of the hard wall model z^ for both meson 
and baryon sectors. 

Now we collect some results from the soft wall model [33] . There were two non-trivial issues 
to be resolved in the original soft wall model. Firstly, so called, the dilaton factor <I> ~ z'^ 
is introduced phenomenologically to explain m^ ~ n. The dilaton factor is supposed to be a 
solution of gravity-dilaton equations of motion. Secondly, the chiral symmetry breaking in the 
model is a bit different from QCD since the chiral condensate is proportional to the quark mass 
in the soft wall model. In QCD, in the chiral limit, where the quark mass is zero, the chiral 
condensate is finite that characterizes spontaneous chiral symmetry breaking. Several attempts 
have made to improve these aspects and to fit experimental values better [65| - |68| . In I66j, a 
quartic term in the potential for the bulk scalar X dual to qq is introduced to the soft wall 
model to incorporate chiral symmetry breaking with independent sources for spontaneous and 
explicit breaking; thereby the chiral condensate remains finite in the chiral limit. Then, the 
authors of [66j parameterized the vev of the bulk scalar Xq such that it satisfies constraints from 
the AdS/CFT at UV and from phenomenology at IR: Xq ~ rUgZ ± az^ as z — )• and Xq ~ z as 



z — )• c«. The constraint at IR is due to the observation 69 that chiral symmetry is not restored 
in the highly excited mesons. Note that Xq ~ z keeps the mass difference between vector and 
axial- vector mesons a constant as z — )• oo. With the parameterized Xq, they obtained a dilaton 
factor $(z) [66]. We list some of results of 66 in Table ^ An extended soft wall model with 



a finite UV cutoff was discussed in 70 , 7l] .In 72 , the authors studied how a dominant tetra 
quark component of the lightest scalar mesons in the soft wall model, where a rather generic 
lower bound on the tetra-quark mass was derived. 
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n 


/9-meson 


p experiment 


ai-meson 


oi experiment 


1 


475 


775.5 


1185 


1230 


2 


1129 


1282 


1591 


1647 


3 


1429 


1465 


1900 


1930 


4 


1674 


1720 


2101 


2096 


5 


1884 


1909 


2279 


2270 


6 


2072 


2149 






7 


2243 


2265 







Table 2: Meson spectroscopy from the modified soft wall model [66]. We show the center values 
of experimental data. In [66^ the experimental data are mostly taken from the particle data 
group |59i, while p(1282) is from |73|. All results are given in units of MeV. 



As long as confinement and non-Abelian chiral symmetry are concerned, the Sakai-Sugimoto 



model 21 , 22 based on a D4/D8/D8 brane configuration (see Appendix O is the only avail 



able stringy model. In this model, properties of light mesons and baryons have been greatly 



studied [2l||22||7§|8T 



In a simple bottom-up model with the Chern-Simons term, it was also shown that baryons 
arise as stable solitons which are the 5D analogs of 4D skyrmions and the properties of the 
baryons are studied [82] , 

5.3 Heavy quarkonium 

The properties of heavy quark system both at zero and at finite temperature have been the 
subject of intense investigation for many years. This is so because, at zero temperature, the 
charmonium spectrum reflects detailed information about confinement and interquark potentials 
in QCD. At finite temperature, due to the small interaction cross section of the charmonium 
in hadronic matter, the charmonium spectrum is expected to carry information about the early 
hot and dense stages of relativistic heavy ion collisions. In addition, the charmonium states 
may remain bound even above the critical temperature Tc- This suggests that analyzing the 
charmonium data from heavy ion collision inevitably requires more detailed information about 
the properties of charmonium states in QGP. Therefore, it is very important to develop a con- 
sistent non- perturb at ive QCD picture for the heavy quark system both below and above the 



phase transition temperature. For a recent review on heavy quarkonium see, for example, 83 

Now we start with the hard wall model to discuss the heavy quarkonium in a bottom-up 
approach. A simple way to deal with the heavy quarkonium in the hard wall model was proposed 
in l84i. Since the typical energy scale involved for light mesons and heavy quarkonia are quite 
different, we may introduce an IR cutoffs z^ for heavy quarkonia in the hard wall model which 
is different from the IR cutoff for light mesons, l/z:^ ~ 300 MeV. Note that in the hard wall 
model there is a one-to-one correspondence between the IR cutoff and the vector meson mass 



1/ Zm ~ rny- In 84 , the lowest vector cc {J/i>) mass, ~ 3 GeV is used as an input to fix the 
IR cutoff for the charmonium, 1/z^ ~ 1.32GeV. With this, the mass of the second resonance is 
predicted to be ~ 7.2 GeV, which is quite different from the experiment m'^ ~ 3.7 GeV. This is 
in a sense generic limitation of the hard wall model whose predicted higher resonances are quite 
different from experiments. Moreover, having two different IR cutoffs in the hard wall model 
may cause a problem when we treat light quark and heavy quark systems at the same time. In 
the soft wall model, the mass spectrum of the vector meson is given by [33] 



ml = 4{n+l)c. (5.9) 

For charmonium system, again the lowest mode (J/ip) is used to fix c, ^/c ~ 1.55 GeV. Then the 
mass of the second resonance ^' is m^i ~ 4.38 GeV, which is 20% away from the experimental 

15 



value of 3.686 GeV [84]. Additionally, the mass of heavy quarkonium such as J/ijj at finite 
temperature is calculated to predict that the mass decreases suddenly at Tc and above Tc it 
increases with temperature. Furthermore, the dissociation temperature is determined to be 
around 494 MeV in the soft wall model J8^ . 

To compare heavy quarkonium properties obtained in a holographic QCD study with lattice 
QCD, the finite-temperature spectral function in the vector channel within the soft wall model 
was explored in [85]. The spectral function is related to the two-point retarded Green function 
by p{u,q) = —2lm.G^{uj,q). The retarded function can be computed following the prescription 
[55] . Thermal spectral functions in a stringy set-up, D3/D7 model, were extensively studied 
To deal with the heavy quarkonium in the soft wall model, two different scales {cp 



m 



86 

and cj/^) are introduced. It is observed in [85] that a peak in the spectral function melts with 
increasing temperature and eventually is flattened at T ~ 1.2Tc. It is also shown numerically 
that the mass shift squared is approximately proportional to the width broadening [85] . Another 
interesting finding in [85] is that the spectral peak diminishes at high momentum, which could 



be interpreted as the J/f/' suppression under the hot wind 87,88 . A generalized soft wall mode 
of charmonium is constructed by considering not only the masses but also the decay constants 
of the charmonium, J/ip and ip' |89J. They calculated the spectral function as well as the 
position of the complex singularities (quasinormal frequencies) of the retarded correlator of the 
charm current at finite temperatures. A predicted dissociation temperature is T « 540 MeV, or 
2.8Tc [89]. 

Alternatively, heavy quarkonium properties can be studied in terms of holographic heavy- 
quark potentials. Since the mass of heavy quarks are much larger than the QCD scale parameter 
^QCD ~ 200 MeV, the non-relativistic Schrodinger equation could be a useful tool to study heavy 
quark bound states. 

^ +V{r)]^{r) = E^{r), (5.10) 



27n. 

where rrir is the reduced mass, rrir = mQ/2. A tricky point with potential models for quarkonia 
is which potential is to be used in the Schrodinger equation: the free energy or the internal 
energy. In the context of the AdS/CFT, there have been a lot of works on holographic heavy 
quark potentials |90| . Hou and Ren calculated the dissociation temperature of heavy quarkonia 
by solving the Schrodinger equation with holographic potentials | 91| . They used two ansatze of 
the potential model: the F-ansatz (U-ansatz) which identifies the potential in the Schrodinger 
equation with the free energy (the internal energy), respectively. With the F-ansatz, J/ip does 
not survive above Tc, while the dissociation temperature of T is (1.3 — 2.1) Tc. For the U- 
ansatz, J/ip dissolves into open charm quarks around (1.2 — 1.7)Tc and T dissociates at about 
(2.5-4.2)Tc. 



We finish this subsection with a summary of the discussion in 20 on the usefulness of Dq/Dp 
systems in studying heavy quark bound states. A Dq/Dp system may be good for ss bound 
states at high temperature since the mesons in the Dq/Dp system are deeply bounded, while 
heavy quarkonia are shallow bound states. However, there exist certain properties of heavy 
quarkonia in the quark-gluon plasma that could be understood in the D4/D6 model such as 
dissociation temperature. 

5.4 Form-factors 

Form factors are a source of information about the internal structure of hadrons such as the 
distribution of charge. We take the pion electromagnetic form factor as an example. Consider 
a pion-electron scattering process vr +6^—7- vr + e~ through photon exchange. The cross 
section of this process measured in experiments is different from that of Mott scattering which 
is for the Coulomb scattering of an electron with a point charge. This deviation is parameterized 
into the pion form factor FT^{q^), where q^ is given by the energy and momentum of the photon 
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UJ 



q . If the pion is a structureless point particle, we have i<V = 1- The pion electromag- 
netic form factor is expressed by, with the use of Lorentz invariance, charge conjugation, and 
electromagnetic gauge invariance, 



where q^ = {p2 — Pi)"^ and J^ is the electromagnetic current, J^ 
charge radius is determined by 



{rl) = 6 



dF^{q^). 



(5.11) 
Y.f(if^nt^'if- The pion 



(5.12) 



In a vector meson dominance model, where the photon interacts with the pion only via vector 
mesons, especially p meson, the pion form factor is given by 



FM^ 



mz 



g2 - impTp{q^ 



m„ 



(5.13) 



Then we obtain the pion charge radius \/{r'^) = vG/nip ~ 0.63 fm. The experimental value 
is ^/{rf) = 0.672 fm 



92 . To evaluate the form factor, we consider the three-point correlation 
function of two axial vector currents which contains nonzero projection onto a one pion state 
and the external electromagnetic current. 



r^a/3(Pl,P2) 



dx / dye'^-'P^-+'^P'^yH0\T{ji,{x)Jpi0)J^,{y)}\0). 



(5.14) 



Alternatively, we can consider two pseudoscalar currents instead of the axial vector currents. The 
three-point correlation function can be decomposed into several independent Lorentz structures. 
Among them we pick up the Lorentz structure corresponding to the pion form factor. 



J^5\P2){P2\JM){P1\J15\0) ^ f^F^{q')p'^p'M+PT. 



(5.15) 



Note that {0\Ja5\p) = ifnPa, where \p) is a one pion state. For more details on the form factor, 
we refer to [93] . 

In a holographic QCD approach, we can easily evaluate the three-point correlation function 
of two axial vector currents (or two pseudoscalar currents) and the external electromagnetic 
current. In \9^, the form factors of vector mesons were calculated in the hard wall model and 
the electric charge radius of the p-meson was evaluated to be (r^) = 0.53 fm^. The number 
from the soft wall model is (r^) = 0.6 
equations predicted (r^) = 0.37 fm^ 



96 



95 



and( 



The approach based on the Dyson-Schwinger 
;) = 0.54 fm2 



97 



The quark mass (or pion 



mass) dependence of the charge radius oTthe p-meson was calculated in lattice QCD: for instance, 

0.55 fm2 



with ttItt 
model 



300 MeV, (H) 



99 



98 



The pion form factor were studied in the hard wall 
and in a model that interpolates between the hard wall and soft wall models [llOOJ. 
The results obtained are \/{r'^) = 0.58 fm [99] and in flOO] yv 



0.58 fm f99] and in flOO] ^/{rl) = 0.500 fm, /(rf) = 0.576 fm, 
depending on their parameter choice. The gravitational form factors of mesons were calculated 
in the hard wall model 



101 102 . The gravitational form factor of the pion is defined by. 



1 



(7r^(p')|e^^(0)K'^(p)) = -5'^\{g^''q^ - q^-q^Qliq^) + AP>'P''Q2{q^)] , 



(5.16) 



where Q^u is the energy momentum tensor, q = p' — p, and P = {p' + p)/2. There are also 
interesting works that studied various form factors in holographic QCD [103f 106 . Form factors 
of vector and axial- vector mesons were calculated in the Sakai-Sugimoto model [107 . 
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Figure 4: QCD phase diagram. 



6 Phases of QCD 

Understanding the QCD phase structure is one of the important problems in modern theoretical 
physics, see |108j for some recent reviews. However, a quantitative calculation of the phase 
diagram from the first principle is extraordinarily difficult. 

Basic order parameters for the QCD phase transitions are the Polyakov loop which charac- 
terizes the deconfinement transition in the limit of infinitely large quark mass and the chiral 
condensate for chiral symmetry in the limit of zero quark mass. The expectation value of the 
Polyakov loop is loosely given by 



(L) 



lim e 



-/3V(r) 



where V{r) is the potential between a static quark-antiquark pair at a distance r, and /3 ~ 
1/r. The expectation value of the Polyakov loop is zero in confined phase, and it is finite 
in deconfined phase. While the chiral condensate, which is the simplest order parameter for 
the chiral symmetry, is non-zero with broken chiral symmetry, vanishing with a restored chiral 
symmetry. Apart from these order parameters, there are thermodynamic quantities that are 
relevant to study the QCD phase transition. The equation of state is one of them. The energy 
density, for instance, has been found to rise rapidly at some critical temperature. This is usually 
interpreted as deconfinement; liberation of many new degrees of freedom. The fluctuations 
of conserved charges such as baryon number or electric charge [109 , llO] is also an important 
signal of the quark-hadron phase transition. The quark (or baryon) number susceptibility, 
which measures the response of QCD to a change of the quark chemical potential is one of such 
fluctuations [l09l[m 



The nature of the chiral transition of QCD depends on the number of quark flavors and the 
value of the quark mass. For pure SU{3) gauge theory with no quarks, it is first order. In the 
case of two massless and one massive quarks, the transition is the second order at zero or small 
quark chemical potentials, and it becomes the first order as we increase the chemical potential. 
The point where the second order transition becomes the first order is called tricritical point. 
With physical quark masses of up, down, and strange, the second order at zero or low chemical 
potential becomes the crossover, and the tricritical point turns into the critical end point. 
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6.1 Confinement/deconfinement transition 

We first discuss the deconfinement transition. In holographic QCD, the confinement to decon- 



finement phase transition is described by the Hawking-Page transition 112 , a phase transition 
between the Schwarzschild-AdS black hole and thermal AdS backgrounds. This identification 
was made in tl2j. One simple reasoning for this identification is from the observation that the 
Polyakov expectation value is zero on the thermal AdS geometry, while it is finite on the AdS 
black hole. See Appendix |F] for some more description of the Hawking-Page transition and the 
Polyakov expectation in thermal AdS and AdS black hole. In low-temperature confined phase, 
thermal AdS, which is nothing but the AdS metric in Euclidean space, dominates the partition 
function, while at high temperature, AdS-black hole geometry does. This was first discovered 
in the finite volume boundary case in |12J. In the bottom-up model, it is shown that the same 
phenomena happen also for infinite boundary volume if there is a finite scale associated with 
the fifth direction [561. 



Here we briefiy summarize the Hawking-Page analysis of 56 done in the hard wall model. 
The Euclidean gravitational action given by 

Sgrav = ~ ^ j ^^^^[^^ Jp.) ' (*^-^) 

where k? = SvrGs and L is the length scale of the AdS^ , there are two solutions for the equations 
of motion derived from the gravitational action. The one is the sliced thermal AdS (tAdS) 

ds"^ = -^ (dr^ + dz"^ + dxl) , (6.2) 

where the radial coordinate runs from the boundary of tAdS space z = to the cut-off Zm- Here 
T is for the compactified Euclidean time-direction with periodicity /3'. The other solution is the 
AdS black hole (AdSBH) with the horizon z^ 

ds' = ^ (/(z)dr2 + -^ + d4^ , (6.3) 

where f{z) = 1 — (z/zh)"^- The Hawking temperature of the black hole solution is T = l/^nzh), 
which is given by regularizing the metric near the horizon. At the boundary z = e the periodicity 
of the time-direction in both backgrounds is the same and so the time periodicity of the tAdS 
is given by 

/3 = 7TZh^/f{ej- (6.4) 

Now we calculate the action density V, which is defined by the action divided by the common 
volume factor of R"^. The regularized action density of the tAdS is given by 

47,3 /•/?' fZiR Ji^ 

and that of the AdSBH is given by 

V2{e) = ^ / dr ^ (6.6) 

1^ Jo Je ^ 



where z = uim{zm, Zh). Then, the difference of the regularized actions is given by 

(6.7) 



( L^-KZh 1 y <r 7^ 






^h 



L^T^Zh f I _ I 



72~ I 7i ^ ) ^m > Zh- 
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When AVg is positive (negative), tAdS (AdSBH) is stable. Thus, at AVg = there exists a 
Hawking-Page transition. In the first case Zm < Zh, there is no Hawking-Page transition and 
the tAdS is always stable. In the second case Zf^ > z^, the Hawking-Page transition occurs at 

n = 2^/V(^^m) , (6.8) 

and at low temperature T < Tc (at high temperature T > Tc) the thermal AdS (the AdS black 
hole) geometry becomes a dominant background. When we fix the IR cutoff by the p meson 
mass, we obtain 1/zm = 323 MeV and Tc = 122 MeV. In the soft wall model, Tc = 191 MeV [5^. 
This work has been extended in various directions. The authors of [ 113| revisited the ther- 
modynamics of the hard wall and soft wall model. They used holographic renormalization to 
compute the finite actions of the relevant supergravity backgrounds and verify the presence 
of a Hawking-page type phase transition. They also showed that the entropy, in the gauge 
theory side, jumps from N^ to A^^ at the transition point 113 . In 114 , the extension was 



done by studying the thermodynamics of AdS black holes with spherical or negative constant 
curvature horizon, dual to a non-supersymmetric Yang-Mills theory on a sphere or hyperboloid 
respectively. They also studied charged AdS black holes [US'] in the grand canonical ensemble, 
corresponding to a Yang-Mills theory at finite chemical potential, and found that there is always 
a gap for the infrared cutoff due to the existence of a minimal horizon for the charged AdS black 



holes with any horizon topology 114 . With an assumption that the gluon condensate melts 



out at finite temperature, a Hawking-Page type transition between the dilaton AdS geometry 



in Eq. (4.4) and the usual AdS black hole has studied in |116|. 

The effect of the number of quark flavors A'^^ and baryon number density on the critical 
temperature was investigated by considering a bulk meson action together with the gravity action 
in |117|. It is shown that the critical temperature decreases with increasing Nf. As the number 
density was raised, the critical temperature begins to drop, but it saturates to a constant value 
even at very large density. This is mostly due to the absence of the back-reaction from number 
density 117 . The back-reaction due to the number density has included in [118 . In |119|, 



deconfinement transition of AdS/QCD with 0{a' ) corrections were investigated. In |120|, 
thermodynamics of the asymptotically-logarithmically-AdS black-hole solutions of 5D dilaton 
gravity with a monotonic dilaton potential are analyzed in great detail, where it is shown that 
in a special case, where the asymptotic geometry in the string frame reduces to flat space with 
a linear dilaton, the phase transition could be second order. The renormalized Polyakov loop in 
the deconfined phase of a pure SU(3) gauge theory was computed in [121 J based on a soft wall 
metric model. The result obtained in this work is in good agreement with the one from lattice 
QCD simulations. 

Due to this Hawking-Page transition, we are not to use the black hole in the confined phase, 
and so we are not to obtain the temperature dependence of any hadronic observables. This is 
consistent with large Nc QCD at leading order. For instance it was shown in [122 123 that 



the Wilson loops, both time-like and space-like, and the chiral condensate are independent 
of the temperature in confining phase to leading order in l/Nc- This means that the chiral 
and deconfinement transitions are first order. The deconfinement and chiral phase transitions 
of an SU(N) gauge theory at large Nc were also discussed in 124[. However, in reality we 



observe temperature dependence of hadronic quantities, and therefore we have to include large 
Nc corrections in holographic QCD in a consistent way. A quick fix-up for this might be to 
use the temperature dependent chiral condensate as an input in a holographic QCD model and 
study how this temperature dependence conveys into other hadronic quantities [125 . 



6.2 Chiral transition 

Now we turn to the chiral transition of QCD based on the chiral condensate. In the hard wall 
model, the chiral symmetry is broken, in a sense, by the IR boundary condition. In case we 
have a well-defined IR boundary condition at the wall z = z^, we could calculate the value of 
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chiral condensate by solving the equation of motion for the bulk scalar X. In the case of the 
AdS black hole we could have a well defined IR boundary condition at the black hole horizon, 
which allows us to calculate the chiral condensate. For instance in |126l|, it is shown that with 
the AdS black hole background the chiral condensate together with the current quark mass is 
zero both in the hard wall and soft wall models. This is easy to see from the solution of Xq in 
the AdS black hole background [126, 127 



Mz) = ^{^rn, ,F,{^-, \, i, ^) + a,z' ^F.i^, ^, | ^)) . (6.9) 

At z = Zfi, both terms in Xq{z) diverges logarithmically, which requires to set both of them zero: 
niq = 0, (T = 0. This is different from real QCD, where current quark mass can be non-zero in 
the regime T > Tc- 



The finite temperature phase structure of the Sakai-Sugimoto model was analyzed in 128| 
to explore deconfinement and chiral symmetry restoration. Depending on a value of the model 
parameter, it is predicted that deconfinement and chiral symmetry restoration happens at the 
same temperature or the presence of a deconfined phase with broken chiral symmetry [128]. 
Phase structure of a stringy D3/D7 model has extensively studied in [129 130 . 



6.3 Equation of state and susceptibility 

Apart from the chiral condensate, various thermodynamic quantities could serve as an indicator 
for a transition from hadron to quark-gluon phase. Energy density, entropy, pressure, and 
susceptibilities are such examples. We first consider energy density and pressure. Schematically, 
based on the ideal gas picture we discuss how the energy density and pressure tell hadronic matter 
to quark-gluon plasma. At low temperature thermodynamics of hadron gas will be dominated 
by pions which are almost massless, while in QGP quarks and gluons are the relevant degrees 
of freedom. Energy density and pressure of massless pions are 

e = 7— ^^ p = ^—T\ (6.10) 

'30 ' -^ '90 ' ^ ' 

where the number of degrees of freedom 7 is three. In the QGP, they are given by 

7r2 7r2 

e = j—T^ + B, p = j—T^-B, (6.11) 

'30 -^ '90 ^ ^ 

where 7 = 37, and B is the bag constant. Apart from the bag constant, the degeneracy factor 
7 changes from 3 to 37, and therefore we can expect that the energy density and pressure will 
increase rapidly at the transition point. Since the dual of the boundary energy- momentum 
tensor T^j^ is the metric, we can obtain the energy density and pressure of a boundary gauge 
theory from the near-boundary behavior of the gravity solution. To demonstrate how-to, we 



follow 131 , 132 . We first rewrite the gravity solution in the Feffer man- Graham coordinate 133 



ds^ = — {gf.^dx^dx" - dz^) . (6.12) 



z 
Next, we expand the metric g^j^y at the boundary 2: — )• 0, 

9^.v = gj^J + z^gj^J + g^l + • • ■ ■ (6-13) 

Now we consider flat 4D metric such that g\i_y = r/^j^. Then g\j,y = and the vacuum expectation 
value of the energy momentum tensor is given by 

(r^,) = const • yg) . (6.14) 
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For example, we consider an AdS black hole in the Feffernian- Graham coordinate 

ds' = ^ [ ^^'ly'^f ' - (1 + ^V4)dx' - dz'^ . (6.15) 

Here the temperature is defined by T = \/^/{'KZh)- Then, we read off 

(T^,) « diag(3/4, 1/4, 1/zl 1/4) , (6.16) 

which satisfies e = 3p. There have been many works on the equations of state for a holographic 
matter at finite temperature 134] . In [135 , the energy density, pressure, and entropy of a 



deconfined pure Yang-Mills matter were evaluated in the improved holographic QCD model [136j. 
The energy density and pressure vanish at low temperature, and at the critical temperature, 
Tc ~ 235 MeV, they jump up to a finite value, showing the first order phase transition. It is 



interesting to note that in [137 some high-precision lattice QCD simulations were performed with 
increasing Nc at finite temperature, and the results were compared with those from holographic 
QCD studies. 

Various susceptibilities are also useful quantities to characterize phases of QCD. For instance. 



the quark number susceptibility has been calculated in holographic QCD in a series of works 126 



138 139 . The quark number susceptibility was originally proposed as a probe of the QCD chiral 



phase transition at zero chemical potential [109, 111 



In terms of the retarded Green function G^j^{io,k), the quark number susceptibility can be 
written as 1140 , 



Xg{T, ^) = - hm ReiG^tico = 0,k)). (6.18) 

In |138| , it is claimed that quark number susceptibility will show a sudden jump at Tc in high 
density regime, and so QCD phase transition in low temperature and high density regime will be 
always first order. Thermodynamics of a charged dilatonic black hole, which is asymptotically 
RN-AdS black hole in the UV and AdS2 xR'^ in the IR, including the quark number susceptibility 



was extensively studied in [141 . The critical end point of the QCD phase diagram was studied 
in [142 by considering the critical exponents of the specific heat, number density, quark number 
susceptibility, and the relation between the number density and chemical potential at finite 
chemical potential and temperature. It is shown that the critical end point located at T = 143 



MeV and n = 783 MeV in the QCD phase diagram 142 



6.4 Dense baryonic matter 

Understanding the properties of dense QCD is of key importance for laboratory physics such as 
heavy ion collision and for our understanding of the physics of stable/unstable nuclei, and of 
various astrophysical objects such as neutron stars. 

To expose an essential physics of dense nuclear matter, we take the Walecka model |143], 
which describes nuclear matter properties rather well, as an example. The simplest version of 
the model contains the nucleon ip, omega meson w, and an isospin singlet, Lorentz scalar meson 
a whose minimal Lagrangian is 

£ = i;{i^ + g^a - g^ ^)i; + -{df^ad^a - mla^) - -F^uPf"" + -mlu^u^ . (6.19) 

Within the mean field approximation, the properties of nuclear matter are mostly determined 
by the scalar mean field a = {ga/rn^)ns and the mean field of the time component of the a; 
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field uq = {g^/'m'^)n, where n is the baryon number density and Ug is the scalar density. For 
instance, the pressure of the nuclear matter described by the Walecka is 



P 



1 

4^ 



p*-ri3 Tr,*2p*^ j_^* 'i^r,(^Z^^PL• 
3 mXr 



(6.20) 



+ J4-^-i4-^ (6-21) 



where 



2 ™* _ ™ .. J^ 

m 



E*F = Jpl + m*^\ m*^ = mN-^n,. (6.22) 



Further, many successful predictions based on the Walecka model and its generalized versions, 
Quantum Hadrodynamics, require large scalar and vector fields in nuclei. This implies that to 
gain a successful description of nuclear matter or nuclei, having both scalar and vector mean 
fields in the model seem crucial. The importance of the interplay between the scalar and vector 
fields can be also seen in the static non-relativistic potential between two nucleons. The nucleon- 
nucleon potential from single cr-exchange and single w-exchange is given by 

V{r) = ^-e-""-"- - 3<Lle-m.r _ ^g_23) 

Att r Att r 

Note that single c-exchange can be replaced by two pion exchange. If ^^^ > g^ and m^ > rria, 



then the potential in Eq. (6.23) captures some essential features of the two nucleon potential 
to form stable nuclear matter: repulsive at short distance and attraction at intermediate and 
long distance. We remark here that the scalar field in the Walecka model may not be the 
scalar associated with a linear realization of usual chiral symmetry breaking in QCD, see for 



instance 1144 



The hard wall model or soft wall model in its original form does not do much in dense matter. 
This is primarily due to its simple structure and chiral symmetry. Suppose we turn on the time 
component of a U(l) bulk vector field dual to a boundary number operator, Vt[z) = fi + pz'^. To 
incorporate this U(l) bulk field into the hard wall mode, we consider U(2) chiral symmetry. The 
covariant derivative with U(l) vector and axial- vector is given by D^X = d^X—iAif^ X+iXAiifj_ 
and it becomes D^X = d^X — iX{AL^ — Ar^). Therefore, the U(l) bulk field V^ = Al^ + tIr^ 
does not couple to the scalar X, meaning that the physical properties of X are not affected by 
the chemical potential or number density. Note, however, that the vacuum energy of the hard 
wall or soft wall model should depend on the chemical potential and number density by the 
AdS/CFT. One simple way to study the physics of dense matter in the hard or soft wall model 
is to work with higher dimensional terms in the action. For instance, the role of dimension six 
terms in the hard wall model was studied in free space [145 . If we turn on the number density 



through the U(l) bulk field, we have a term like X^Fy, where Fy is the field strength of the 



bulk U(l) gauge field 146 . Then we may see interplay between number density and chiral 



condensate encoded in Xq. In 147 , based on the hard wall model with the Chern-Simons term 
it is shown that there exists a Chern-Simons coupling between vector and axial-vector mesons 
at finite baryon density. This mixes transverse p and oi mesons and leads to the condensation 
of the vector and axial-vector mesons. The role of the scalar density or the scalar field in the 
hard wall model was explored in |148| . In 149| , a back-reaction due to the density is studied in 



the hard wall model. 

Physics of dense matter in Sakai-Sugimoto model has been developed with/without the 
source term for baryon charge [150] . For instance, in [151 localized and smeared source terms 



are introduced and a Fermi sea has been observed, though there are no explicit fermionic modes 
in the model. A deficit with the Sakai-Sugimoto model for nuclear matter might be the absence 
of the scalar field which is quite important together with U(l) vector field. The phase structure 
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of the D3/D7 model at finite density is studied in [152 . The nucleon-nucleon potential is playing 



very important role in understanding the properties of nuclear matter. For example, one of the 
conventional methods to study nuclear matter is to work with the independent-pair approxima- 
tion, Brueckner's theory, where two-nucleon potentials are essential inputs. Holographic nuclear 



forces were studied in 153-156 



7 Closing remarks 

The holographic QCD model has proven to be a successful and promising analytic tool to study 
non-perturbative nature of low energy QCD. However, its success should always come with 
"qualitative" since it is capturing only large Nc leading physics. To have any transitions from 
"qualitative" to "quantitative", we have to invent a way to calculate subleading corrections in 
a consistent manner. A bit biased, but the most serious defect of the approach based on the 
gauge/gravity duality might be that it offers inherently macroscopic descriptions of a physi- 
cal system. For instance, we may understand the QCD confinemnt/deconfinement transition 
through the Hawking-Page transition, qualitatively. Even though we accept generously the 
word "qualitatively" , we are not to be satisfied completely since we don't know how gluons and 
quarks bound together to form a color singlet hadron or how hadrons dissolve themselves into 
quark and gluon degrees of freedom. In this sense, the holographic QCD can not be stand-alone. 
Therefore, the holographic QCD should go together with conventional QCD-based models or 
theories to guide them qualitatively and to gain microscopic pictures revealed by the conven- 
tional approaches. 

Finally, we collect some interesting works done in bottom-up models that are not yet properly 
discussed in this review. Due to our limited knowledge, we could not list all of the interesting 
works and most results from top-down models will not be quoted. To excuse this defect we refer 



to recent review articles on holographic QCD 1157 



Deep inelastic scattering has been studied in gauge/gravity duality [158 -169 . Light and 
heavy mesons were studied in the soft- wall holographic approach flTO]. 

Unusual bound states of quarks are also interesting subjects to work in holographic QCD. 
In [171 , the multi-quark potential was calculated and tetra-quarks were discussed in AdS/QCD. 



Based on holographic quark-antiquark potential in the static limit, the masses of the states 
X(3872) or Y(3940) were predicted and also tetra-quark masses with open charm and strangeness 
were computed in [172] . A hybrid exotic meson, 7ri(1400), in [173]. The spectrum of baryons 



with two heavy quarks was predicted in 174[. 



Low-energy theorems of QCD and spectral density of the Dirac operator were studied in the 
soft wall model [175[. A holographic model of hadronization was suggested in 176 . 



The equation of state for a cold quark matter was calculated in the soft wall metric model 



with a U(l) gauge field. The result is in agreement with phenomenology 177 . 
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A Bulk mass and the conformal dimension of boundary opera- 
tor 

In this Appendix, we summarize the relation between the conformal dimension of a boundary 
operator and the bulk mass of dual bulk field. We work in the Euclidean version of AdSd+i, 

d'^" = A-2Y.(d^n'- (A.l) 

A.l Massive scalar case 

We first consider a free massive scalar field whose action is given by 



S = ^fd''+^x^{df,<pd''^ + m^(f>^). (A.2) 



Let the propagator of 4> be K{x^,x;x'). To solve for cf) in terms of its boundary function (f>Q, 
we look for a propagator of 4>, a solution K{x^,x;x') of the Laplace equation on B^+i whose 
boundary value is a delta function at a point P on the boundary. We take P to be the point at 
xo — )• oo. The boundary conditions and metric are invariant under translations of the Xi, then 
we can consider K as a function of only xq, thus K{x^,x; P) = K{x^). Then, the equation of 
motion is 

-(-°)^^^^(-°)-^^^^ +-^) K{x^) = 0, (A.3) 

where we used 

-La,V58'- = (.»)" ;^(.°)-«^. (A.4) 

We analyze the equation of motion near the boundary, x'' — )• 0, and take K[x^) oc (x^)^^*^. 
From the equation of motion, we have 

-(A + d)A + m2 = 0, (A.5) 

where A is the larger root A = A+. The conformal dimension A of the boundary operator is 
related to the mass m on AdS^+i space by A = d + A+. Thus, we obtain 

(A - d)A = m^ , (A.6) 

or 

A = ^ (d + ^d^ + 4m2) . (A.7) 

A.2 Massive p-form field case 



Consider a massive p-form potential 178 



A=\A^^...y,^dx^^---dx^^ . {k.l 



The free action of A is 



S=\l {TA*T+ Tr?A A M) , 

2 JAdSa+i 

where J- = dA is the field strength p+ 1 form. The variation of this action is 
dS= [ {-{-l)P6AAd*\F + m'^6AA*A) , 



(A.9) 
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and then the classical equation of motion for A from (A. 9) is 

{-l)Pd*dA-m^*A = 0. 



(A.IO) 



In addition, A satisfies d*A = 0. By using the metric (A.l), the equation of motion (A.IO) can 
be written as 



[{x^dl -{d+1- 2p)x°9o + {d+l-2p- m^)] Aoi,..A, 
[{xyOl -{d-1- 2p)x^do - m^] A,..., 

= 23;° {di^uJoi2...ip + {-lY~'^di^ujoi3...iph H ) • 

Now from the vielbein Ca = x^6a, we introduce fields with flat indices 



0, 



Then the equations of motion (A. 11) of Aoi2...ip becomes 

[{xydl -id- l)x05o - (m2 + p2 _ pd)] Ao^,..,p = 0. 
We consider 



A 



012 ■■■ip 



(x") 



0\-A 



as X —7- 0. Then substituting this in (A. 14) gives 



= [{x^fdi -{d- l)x°5o - (m^ + p2 _ pd)] (x^)-^ 

= "(x°)29o(-A(x°)-^-i) -{d- l)a;0(-A(x°)-^-i) - (m^ + / _ pd){x^) 

= [A(A + 1) + A(d - 1) - (m^ + / _ pd)] (x^)-^ 

= [X{X + d)-{m'^+p'^-pd)]{x^)-^ 

and therefore we obtain the relation 

X{X + d) =rr? +p^ -pd. 



0\-A 



With A = d + A, we have 



and we finally arrive at 



or 



(A - (i) A = m^ +p^ -pd 
^ (A - p)p + (A - p)(A - d) = m^ 

(A -d + p)(A-p) = m^, 
A = Vd + v/('i - 2p)2 + 4.rrA . 



A. 3 General cases 



(A.ll) 
(A.12) 



(A.13) 

(A.14) 
(A.15) 



(A.16) 



(A.17) 
(A.18) 



Now for completeness, we list the relations between the conformal dimension A and the mass 
for the various bulk fields in AdSd+i- 



1. scalars |3j: A± = \{d ± yj d^ + Am?) , 

2. spinors ^: A = ^(d + 2|m|), 

3. vector^ A± = \{d ± ^ {d - 2)"^ + Am?) , 
^Entries 3. and 4. are for forms with Maxwell type actions. 
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4. p-forms 



178 



: A± = i(d± v^(d-2p)2 + 4m2), 

5. first-order ((i/2)-forms {d evenfj A = ^{d + 2|m|), 

6. spin-3/2 fTsTl [l82|: A= ^(d + 2|m|), 



7. massless spin-2 fISS] : A = d. 

B D3/D7 model and U(l) axial symmetry 

In the original AdS/CFT, the duahty between type IIB superstring theory on AdS^ x S^ and 
A/" = 4 super Yang-Mills theory with gauge group SU{Nc) can be embodied by the low-energy 
dynamics of a stack of Nc D3 branes in Minkowski space. All matter fields in the gauge theory 
produced by the D3 branes are in the adjoint representation of the gauge group. To introduce 
the quark degrees of freedom in the fundamental representation, we introduce some other branes 
in this supersymmetry theory on top of the D3 branes. 

B.l Adding flavour 

It was shown in 



17 



that by introducing A'^^ D7 branes into AdS^ x S^, Nj dynamical quarks can 
be added to the gauge theory, breaking the supersymmetry to AA = 2. The simplest way to treat 
D3/D7 system is to work in the limit where the D7 is a probe brane, which means that only a 
small number of D7 branes are added, while the number of D3 branes Nc goes to infinity. In 
this limit Nf <C Nc we may neglect the back-reaction of the D7 branes on AdS^ x S^ geometry. 
In field theory side, this corresponds to ignoring the quark loops, quenching the gauge theory. 
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Table 3: The D3/D7-brane intersection in 9 + 1 dimensional flat space. 



The D7 branes are added in such a way that they extend parallel in Minkowski space and 
extend in spacetime as given in Table [3} The massless modes of open strings both ends on the 
Nc D3 branes give rise to AA = 4 degrees of freedom of supergravity on AdS^ x S'^ consisting 
of the SU{Nc) vector bosons, four fermions and six scalars. In the limit of large Nc at fixed 
but large 't Hooft coupling A = Qym^c = QsNc ^ 1 , the D3 branes can be replace with near 
horizon geometry that is given by 



ds^ 



m 
m 



(-dt^ + dxj + dxj + dxl) + 



— . 



{-df + dxl + dxl + dxl) + -^ {dp^ + P^d^l + dyl + dy|) 



(B.l) 



R is the radius of curvature 



where y = (yi, . . . , ye) parameterize the 456789 space and r^ 

R^ = y/AirggNca' and d^l is the three-sphere metric. The dynamics of the probe D7 brane is 



f. 



described by the combined DEI and Chern-Simons actions 5,184 , 



Sd7 = -T7 f d^x V-det {P[gU + ^ira'Fab) + ^^^^Ty f P[C^^^] A F A F . (B.2) 



where g is the bulk metric (B.l) and C*-^' is the four form potential. T7 = l/((27r)^g(sa'4) is the 



D7 brane tension and P denotes the pullback. Fab is the world-volume field strength. 



See 180 for d — A case. 
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The addition of D7 branes to this system as in Table [3] breaks the supersymmetry to J\f = 2. 
The hghtest modes of the 3-7 and 7-3 open strings corresponds the quark supermultiplets in 
the field theory. If the D7 brane and the D3 brane overlap then S0(6) symmetry is broken into 
5*0(4) X SO(2) ~ SO{2)pi x S0{2)l x C/(1)_r in the transverse directions to D3 and so preserves 
1/4 of the supersymmetry. The S0{4) rotates in 4567, while the 50(2) group acts on 89 in|3J 
The induced metric on D7 takes the form, in general, as 

dslj = ^r^^udx^'dx^ + ^^ ((1 + v'i + viW + p'dVif) (B.3) 

where 2/5 = dy^/dp and y'^ = dye/dp. When the D7 brane and the D3 brane overlap, the 
embedding is 

ys = 0, 2/6 = (B.4) 

and the induced metric on the D7 brane is replaced by 

2 r>2 

dslj = j^Vf.udx^dx'' + -J {dp^ + p^dnl) . (B.5) 

The D7 brane fills AdS^ and is wrapping a three sphere of S^. In this case the quarks are 
massless and the R-symmetry of the theory is SU{2)ji x U{1)r and we have an extra C/(l)/j 
chiral symmetry. 

If the D7 brane is separated from the D3 branes in the 89-plane direction by distance L, then 
the minimum length string has non-zero energy and the quark gains a finite mass, niq = L/2TTa' . 
It is known that the R-symmetry is then only SU{2)ji and separation of D7 and D3 breaks the 
50(2) ~ U{1)r that acts on the 89-plane. In this case, we can set for the embedding as 

ys = 0, ye = yeip) ■ (B.6) 

Then, the action for a static D7 embedding (with Fab zero on its world volume) becomes 

Sd7 = -T7 j d^xV-det {P[gU) = -T7 j d^xyj-detgabyjl + g^^dayidby^gij 

= -Tj I d^xes^y 1 + {dpy^y + {dpy^Y , (B.7) 

where i, j = 5,6 and £3 is the determinant from the three sphere. The ground state configuration 
of the D7 brane is given by the equation of motion with 2/5 = 



dp 



P^dpye 



0. (B.^ 



The solution of this equation has an asymptotic behavior at UV (p — )• 00) as 

2/6 ^ m + 4 + • • • • (B.9) 

p^ 



Now we can identify [158 that m corresponds to the quark mass and c is for the quark condensate 
{ipil^) in agreement with the AdS/CFT dictionary. 

B.2 Chiral symmetry breaking 

One of the significant feature of QCD is chiral symmetry breaking by a quark condensate '^p^lJ. 
The C/(l) symmetry under which ip and ip transform as '(/' — ?• e~^°^ip, ip — >• e*°^ in the gauge 
theory corresponds to a U{1) isometry in the 2/52/6 plane transverse to the D7 brane. This U{1) 
symmetry can be explicitly broken by a non- vanishing quark mass due to the separation of the 

28 



D7 brane from the stack of D3 branes in the y^ + ii/q direction. Assume the embedding as 2/5 = 
and ye ~ c/p^, then by a smah rotation e~*^ on 2/5 + iyQ generates 2/5 ~ ec/p^ and yg ~ yQ up to 
the ©(e^) order. 

In [Is] the embedding of a D7 probe brane is embodied in the Constable-Myers background 
and the regular solution yg ~ "^ + c/p^ of the embedding y^ = 0, ye = Veip) shows the behavior 
c 7^ as m —7- which corresponds to the spontaneous chiral symmetry breaking by a quark 
condensate. In |l85| , the chiral symmetry breaking comes from a cosmological constant with a 
constant dilaton configuration which is dual to the AA = 4 gauge theory in a four-dimensional 
AdS space. 

B.3 Meson mass spectrum 

The open string modes with both ends on the flavour D7 branes are in the adjoint of the U{Nf) 
flavour symmetry of the quarks and hence can be interpreted as the mesonic degrees of freedom. 
As an example, we discuss the fluctuation modes for the scalar fields (with spin 0) following the 
argument of [186 . The directions transverse to the D7 branes are chosen to be ys and ye and 



the embedding is 

y5 = + x, y6 = i^ + '^ (B.IO) 

where Sy^ = x and 5ye = (p are the scalar fluctuations of the transverse direction. To calculate 
the spectra of the world- volume fields it is sufficient to work to quadratic order. For the scalars, 
we can write the relevant Lagrangian density as 



£d7 = -T7^/-detP[g] 



ab 



T7v^-detyaf,-y/l + g''^ {daXdbXabb + da^db^9&&) 



02 

-Tyv^-detsfabA/l + a^'^^T {daXdbX + da^pdbf) 






~ -Try/-detgab U + —g"' {daXdbX + da^db'p)] (B.ll) 

where P[g]ab is the induced metric on the D7 world-volume. In spherical coordinates with 
r^ = p'^ + L'^, this can be written as 

Cm ^ -Trp^e^ (l + ^^^g-b {daxdbX + da^dk^)\ (B.12) 

where £3 is the determinant of the metric on the three sphere. Then the equations of motion 
becomes 

9a(^^<7"''5.$)=0, (B.13) 

where $ is used to denote the real fluctuation either x or '■P- Evaluating a bit more, we have 

_^_d,d^^ + ^^P(P'^P^) + ^V.V*ci> = , (B.14) 

where Vj is the covariant derivative on the three-sphere. We apply the separation of variables 
to write the modes as 

$ = (j){p)e^^-'^y\S^) , (B.15) 

where y^{S'^) are the scalar spherical harmonics on S^, which transform in the {£/2,i/2) repre- 
sentation of SO (4) and satisfy 

v*v»y = -i{i + 2)y . (B.16) 

The meson mass is defined by 

M^ = -k'^. (B.17) 
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Now we define q = p/ L and M^ = —k'^R^/L'^, and then the equation for (pip) is 



dj,cb + -d,6 + 



M' 



i{£ + 2) 



(1 + ^^ 



i2^2 



. 



(B.18) 



This equation was solved in 186 in terms of the hypergeometric function. To solve the equation, 
we first set 

HQ) = 0'{^ + QY"Pie), (B.19) 

where 

2a = -l + \/l + M2 > 0. (B.20) 



With a new variable y = —g, (B.18) becomes 



y(l - y)P"{y) + [c-{a + b + l)y]P'iy) - ahP{y) = 



(B.21) 



where a = —a, b = —a + i + 1 and c = i + 2. The general solution is taken by a > and by 
noting that the scalar fluctuations are real for — oo < y < 0, one finds, up to a normalization 
constant, the solution of d) 



^ip) 



P 



-a, 



-a + £+l; i + 2; -p^/L^). 



{p^ + L^Y 
Imposing the nomalizability at p — )• oo, we obtain 

-a + £+l = -n, n = 0,1,2,.... 

The solution is then 



<^{P) 



(p2 + L2) 



n+i+ 



J F (-(n + ^ + 1) , -n ; £ + 2 ; -p^/L^ 



and from the condition (B.23) we get 

'^ -A{n + i+l){n + l + 2) . 



M^ 



(B.22) 



(B.23) 



(B.24) 



(B.25) 



Then by the definition of meson mass (B.17), we derive the four-dimensional mass spectrum of 
the scalar meson 



2L I 

^sin.i) = j^y^{n + e+l)in + i + 2) . 



(B.26) 



B.4 Mesons at finite temperature 



In previous sections, we have focused on gauge theories and their gravity dual at zero tem- 
perature. To understand the thermal properties of gauge theories using the holography, we 
work with the AdS-Schwarzschild black hole which is dual to AA = 4 gauge theory is at finite 
temperature [3,12 . The Euclidean AdS-Schwarzschild solution is given by 



ds' 



K{r) 



dT^ + W 



dr"^ 



K{r) i?2 



+ -2 df2 + R^dOi , 



where 



K{r) 



r^ I 1 



4 
»,4 



(B.27) 



(B.28) 



For r ^ rff, this approaches AdS^ x S^ and the AdS radius R is related to the 't Hooft coupling 
by R^ = v4vrAa'. Note that the S^ parameterized by r collapses at r = rn, which is responsible 
for the existence of an area law of the Wilson loop and a mass gap in the dual field theory. 
This geometry is smooth and complete if the imaginary time r is periodic with the period 



30 



/? = E?"K/rH- The temperature of the field theory corresponds to the Hawking temperature is 
given by T = 1//3 = r///(i?^7r). At finite temperature, the fermions have anti-periodic boundary 
conditions in r direction |12| and the supersymmetry is broken. In addition, the adjoint scalars 
also become massive at one loop. Thus, fermions and scalars decouple. 

We now introduce D7 branes in this background. It is convenient to change the variable in 
the metric (B.27) such that the it possesses an explicit flat 6-plane. To this end, we change the 



variable from r to if as 



r dr 



K{r) 



' H 



dw 



w^ 



and take R= 1. One of the solutions is 



2w' 



r' + 



' H 



or r 



4 I 4 
W + Wff 



(B.29) 



(B.30) 



and wh = ^hI^I^- Then the metric (B.27) becomes 



ds' 



w'^{w'^ + wjj 



"'"^' dt^ + (w' + %) dx2 + -^{dw-^ + w'dni) 



''HJ 



{w^ 



'H) 



w'^{w'^ + wjj 



^" dt^ +{w^ + ^\ dx^ + -^{dp'' + p^dn'i + dw'i + dw'^) 



(B.31) 



where df]| is the three-sphere metric. The D7 brane is embedded in the static gauge, which 
worldvolume coordinates are now identified with xq, 1,2,3 and 11^1,2,3,4 and the transverse fluctua- 
tions will be parameterized by 1^5 and wq. The radial coordinate is given by 



w 



i=l 



2 2 I 2 I 2 

Wi = P +W^+Wq. 



(B.32) 



At large w, this geometry asymptotically approaches AdS^ x S^ and the D7 embedding should 
approach the constant solutions W5 = 0,wq = const which is the same as the exact solution in 
section B.l To consider the deformation, we take the following ansatz for the embedding 

■^5 = 0, ■wq = wq{p) (B.33) 

Then, the action of D7 becomes 

Sd7 = -Tj / d^xyJ-d.(iiP[g\ab = -Tj / d^xyj -deigab\jl + 9°-^ {daW^dbW^g^^ + daWQdbUiQgQe) 

= -T7 f d^X €3 gip,W5,W6)^l + {dpW5y + (OpWe)^ , (B.34) 

where €3 is the determinant of the three-sphere metric and the function g{p, W5,wq) is given by 



g{p,W5,WG) 



(u;4 



0^(^^ + 0V_^3,i 



,16 



W 



H 



With the assumption ^5 = 0, the equation of motion takes the form 



d^ 
dp 



pU 1 



w 



H 



dpWe 



(p' + wlYJ VI + {dpWef 



8wj^p wq 



l + idpW6y = 0. 



(B.35) 



(B.36) 



With the solution of this equation, the induced metric on the D7 brane is given by 



ds 



D7 



{w 



rA 



w 



%? 



w-^{w^ + W]^^ 



w 



df + ( u,^ + ^^ dx^ + i±i?f^d,^ + ^,dnl , 



w^ 



w 



(B.37) 



31 



with w'^ = p^ + 'Wq{p) and the D7-brane metric becomes AdS^ x S"^ for p ^ 
The asymptotic solution at large p is of the form 



Wj^,Wq. 



Wq{p) r^ 101 + 



P- 



,2 • 



(B.38) 



As discussed in section B.l , the parameters m and c are interpreted as a quark mass and bilinear 
quark condensate {ip^), respectively. With suitable boundary conditions for the second order 



equation ( |B.36 ), one can solve it numerically by the shooting method. The results are shown in 
Figure [5j 



We(p) 




Minkowski 
embeddings 
( m > iric ) 

■^^m = iric 

Black hoie 
embeddings 
( m < mc) 



Figure 5: Two classes of regular solutions in the AdS black hole background. 



We see that there can be two different classifications for the D7 brane embeddings. First, 
for large quark masses the D7 brane ends outside the horizon. It can be interpreted that the D7 
brane tension is stronger than the attractive force of the black hole. Such a D7 brane solution 
is called a Minkowski embedding. They behave similarly to the supersymmetric solutions in 
AdS^ X S^. Second, for small masses the D7 brane ends at the horizon w = wh at which the S^ 
of the black hole geometry collapses. This is called a black hole embedding. These two classes of 
embeddings also differ by their topology. The D7 brane topology is i?^ x i?^ x S^ for Minkowski 
embedding and R^ x S'^ x B^ for black hole embedding. The change in the topology is related 
to a phase transition in the dual field theory. 

C D4/D8/D8 model and non-Abelian chiral symmetry 



The D3/D7 system is a supersymmetric configuration which gives gauge theories in the ultra- 
violet and only has a U{1)a symmetry. To realize more realistic non-Abelian chiral symmetry, 
Sakai and Sugimoto proposed D4/D8/D8 brane configuration 21,22 which is intrinsically non- 
supersymmetric. They placed probe D8 and D8 branes into the Nc D4 background of the 
fundamental representation of the SU{Nc) gauge group. The D4/D8/D8 intersection in 9+1 
dimensional flat space is given in Table |4j 
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o 



Table 4: The D4/D8/D8 brane intersection in 9 + 1 dimensional flat space. 
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In this configuration, Nc D4-branes are compactified on a direction (x^) wrapped on S^. 
We impose anti-periodic boundary conditions on this S^ in order to break the supersymmetry. 
Nf D8 and D8 pairs are put in the transverse to the S^. Then, from the D4-D8 and D4-D8 
open strings we obtain Nj flavors of massless chiral and anti-chiral quark fields of the U{Nc) 
gauge group. If we take the strong coupling limit to take the large Nc dual, the D8 and D8 
will overlap into a single curved D8 brane, which is interpreted as non-Abelian chiral symmetry 
breaking. The U{Nf)Ds x U{Nf)j^ gauge symmetry of the D8 and D8 branes is interpreted as 
the U{Nf)i X U{Nf)ji chiral symmetry. These configurations are sketched in Figure ml 




D4 




Figure 6: Sketch of the D4/D8/D8 configuration. 



C.l Background D4 and probe D8 branes 

In order to obtain a holographic dual of the the large Nc gauge theory with non-Abelian chiral 
symmetry, we consider the SUGRA description of the D4/D8/D8 system as discussed above. 
Assuming Nf <^ Nc, we treat D8-D8 pairs as probe D8 branes embedded in the D4 background. 
By taking the near horizon limit of the geometry of Nc stack of D4 branes wrapped on a circle, 
we obtain the D4 brane solution geometry 

ds' = (I) '^' {v^.dx^dx'^ + f{U)dr') + {§y\^^+ U'dni^ (C.l) 

with f{U) = 1 — Uf^j^/U^. Here x^ (/i = 0,1,2,3) and r are the directions along which the 
D4-brane is extended and R is the radius of curvature related to the string coupling gg and 
string length Is as R"^ = irggNcls- There is a nonzero four- form flux F4 = dC^ = 2TTNce4/V4 
with the volume form €4 and a dilaton e"*^ = gs {U / R)~ . U is the holographic direction and 
bounded from below by the condition U > Ukk- As [/ — )■ Ukk the radius of S^ parameterized 
by T shrinks into zero and then the D8 and D8 branes are connected at some point U = Uq. 
In this case the gauge symmetry is U{Nf). The configuration is then interpreted as a theory 
with massless quarks with chiral symmetry breaking, which is sketched on the right hand side 
of Figure [6} In order to avoid a singularity at U = Ukk-, t must be periodic with period 
St == 47ri?3/2/3[7Vj^ ^j^j^j^ defines the Kaluza-Klein mass, Mkk = 2-k/5t = 3C/J/J/2i?3/2. 
Next, we consider the induced metric on the D8 probe brane in the D4 background with an 
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ansatz U = U{t). Then, dll = [dU /dr) t = U'{t) dr and the induced metric is given by 



UV'^ ^ ^^ , (UV'^ ^,^^, (rV'^ U'^\,2 fRV^' 



dsU = [^J^j V^udx^^dx^ + ( (^- j f{U) + (^- j -^1 dr' + (^- j U'd% , (C.2) 
The D8 brane DBI action now becomes 



Sds ^ j d^xdre^e-^^J-detigns) ^ j d^ xdrU* J f {U) + (^\ j^ (C.3) 



Since the integrand of (C.3) does not exphcitly depend on r, we obtain the energy conservation 



" ' '',w -c] = 4- („„, "'EL,.„, ] = 0. (C.4) 



dr \dU' J dT\f{U)+K^U''^/U'^f{U) 



Let ?7(0) = C/q and assume t/'(0) = at r = 0, then the solution of (C.4) wih be 



r{U) = Ulf{U^f'^ r "^^ ^ =. (C.5) 

Juo{U/Rf/^f{U)^JWjm^WW) 

It can be shown that r(C/ — >• oo) is a monotonically decreasing function of Uq, roughly r(oo) ~ 

Uq , that is varying from t{oo)\u=Ukk ~ ^/^ ^o t{oo)\u^oo = 0. In the limit Uq = Ukk, the 
D8 and D8 branes are at antipodal points on the S^ parameterized by r. And when Uq — )■ oo, 
the D8 and D8 separate far away from each other. Here we concentrate on the case Uq = Ukk- 
For the sake of convenience, we perform the change of variables 

U^ = Uf,K + UKKr\ 0=f^r='^r. (C.6) 



Now, the induced metric (C.2) on a D8 brane becomes 



dsls =[lf V,.dxV + ^ (I) '^' (^rfr^ + r'de^ + (f ) '^' U'd^j 



(C.7) 



We change the variables once more as 

y = rcos9, z = rsm9, {CI 



accordingly, the dr and d9 part of (C.7) read 



Ukk J 2 , „2,/)_/i ^,/„^^2^T^2 , n i./„\„.2\,.2 



u 



-dr' + r'de = (1 - h{r)z')dz' + (1 - h{r)y')dy' - 2h{r)zydzdy (C.9) 



with h{r) = l/r^(l — Ukk/U). Near U = Ukk, (C.7) approaches a flat two-dimensional plane 
and h{r) is a regular function in the neighborhood of r = 0. Then, y{x'^, z) = can be a solution 
of the equation of motion of the probe D8 brane's world- volume theory and the stability of the 
solution can be shown by examination of small fluctuations around it f21]. Thus, from (C.6) we 
have 

U^ = UI,k + Ukkz^ (CIO) 

and the induced metric on the D8 brane 
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C.2 Gauge field and meson spectrum 

Now, we consider the gauge field on the probe D8 brane configuration. The gauge field on the D8 
brane Am {M = 0, 1, 2, 3, 5, 6, 7, 8, z) has nine components, and among them we are interested 
in the 50(5) singlet states, so we can set Aa = for a = 5, 6, 7, 8. We assume that A^ and Az 
are independent of the coordinates on the S . Then the DBI action for D8 branes becomes 



Sd8 = -Tg / (fxy^-det {gab + 2-Ka'Fab) + Scs 



T{2TTa'Y / (Txdz 



AU 



o UkK 



+ 0{F^), (C.12) 



where f = {2/?,)R^''^U jl^T^VA^Qg ^. We expand the gauge fields A^ and A^ in terms of the com- 
plete sets of the Kaluza-Klein (KK) modes profile functions {'(/'^(-z)} and {(f)n{z)}, respectively, 



n>0 

Az{x,z) = Y,ip^''\x)<Pn{z). 



n>0 



and consequently, the field strengths read 

F,,{x,z) = Y, {d^v^:'\x)Mz)) = j;FW(:r)^„(, 

n>0 n>0 



(C.13) 
(C.14) 

(C.15) 
(C.16) 



n>0 



Then, the action (C.12) becomes 



Sds = — T{2TTa) I d xdz y 



n,m>0 



AU 



^^n)i.MKn)^^^^ 



For convenience, we introduce a new variable Z and the function K{Z) as 

Z = z/Ukk , and K{Z) = 1 + Z^ 



(C.17) 



(C.18) 



With these, (C.17) changes into 



5^8 = - f{2T:a'fR^ / d'^xdZ 



n,m>0 



^i^-i/3(^)^Wi7/^Kn)^^^^ + iM|:^i^(Z)t;Wt;'^(")az^n9zV'. 



T{2Tra'fR^ / d^x dZ -M]^j^K{Z) 



E {uKK^^.^^''^^^'^'^'^^n(t>^ - 2C/|.;^a^v^(")i;^('")0„5z^« 



n,r?i>0 



Now we choose ipni^) (n > 1) as the eigenfunctions 



(C.19) 



(C.20) 



35 



with the normalization condition 



f{2Tra')^R^ f dZR-^/^tl^m^n = Sn 



to arrive at 



f{27ra'fR^ j dZKdz^mdzi^n = K^n 
Then, the action (C.19[) becomes 



'' n>0 



.l^H^MKn) _ ^Mi^A„T;Wz;'^(") 



-,{n) 



+ if"-"-' parts 



(C.21) 

(C.22) 

(C.23) 



Similarly, we can normalize the profile fmictions {4>n\ of ip^^'. By observation (C.22), we can 
choose 

and (j)o = -T^. (C.24) 



1 



K = —rj^ dzipn (ra > 1) 



K{zy 



for some constant C. Then from (C.22) the orthonormal condition for {(pn} becomes 

f(2W)2M|.^i?3 [ dZK<i)^4>n = 5^mn. (C.25) 



The definitions in (C.24) imply that 0o is orthogonal to (/>„ for n > 1 since j dZ K{Z)(j)Q(f)n ~ 
f dZ dzi/^n = from (C.25). Taking n = m = in (C.25), we determine C as C~^ = 
f{2T^a'fMJ^j^R^TT since J dZ K{Zy^ = f dZ {1 + Z^)"! = vr. Note that Vo can be deter- 
mined as 

V'o(^) = / dz-^ = Ctan'^z (C.26) 

Jo 1 + ^ 

which is not normalizable. However, the field strength is normalizable and ipo can be considered 

as the zero mode of the eigenfunction equation (C.20). Now from (C.24), F^z becomes 



F^z{x, z) = Y, [d^^^'^Un - <)9,V« 

n>0 

n>l 

We can absorb An M^j^d^tp^^' into t'/T , and then we obtain 



(C.27) 



(C.28) 



n>l 



Then again with the change of variable (C.18), the action (C.12) can be written as 
5^8 = - f{2-Ka'fR^ [ d^x dZ 

- E 

n,m>0 



^i^-i/3(^)^Wi7/^Kn)^^^^ + lMl,j,K{Z)v^r^\^^'^)d^^^d^^^ 



f{2TTa'fR^ I d^x dZ ^Mi^K(Z)?7i^a^v9(0)9'^c^(0)(/)2. 



(C.29) 



With the normalization conditions (C.21), (C.22), and C — )■ C/Ukk, the action (C.29) becomes 



5*08 = d X 



-5^(^(0) 9^V9(°) + Y, f-^^/."^^'^'^^"^ - -mSz;(")t;A'(-) 



n>l 



r)"''n"ll 



(C.30) 
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where 

ml = MIkK . (C.31) 

The KK modes of A^ are t"^" {n>l) and they are regarded as the massive vector meson field. 
Also, we interpret ^f^^', which is the KK mode of A^, as the pion field or the Nambu-Goldstone 
boson of the chiral symmetry breaking. 



By solving (C.20) numerically, we find the eigenvalues corresponding to the masses of vector 
mesons 

A„ = 0.67, 1.6, 2.9, 4.5, ..., (C.32) 

and we compare the meson mass ratio obtained in this model with the experimental data 



2Ai — y -^ = 2.51, 

2 



Ai ■ m'p 



^' - 4.3 ^ ^^^ = 3.56. (C.33) 



A 



D Surface gravity of a Schwarzschild black hole 

The surface gravity is the gravitational acceleration experienced by a test body (with negligible 
mass) close to the surface of an object. For a black hole, the surface gravity is defined as 
the acceleration of gravity at the horizon. The acceleration of a test body at a black hole 
event horizon is infinite in relativity, therefore one defines the surface gravity in a different way, 
corresponding to the Newtonian surface gravity in the non-relativistic limit. Thus for a black 
hole, the surface gravity is defined in terms of the Killing vector which is orthogonal to the 
horizon and here its event horizon is a Killing horizon. For the Schwarzschild case this value is 
well defined. Alternatively, one can derive the same value as a period of the imaginary time in 
the Euclidean signature, we will discuss both points of view. 

D.l Orthogonal Killing Vector 

The horizon of a black hole is a null surface. It means that any vector normal to the surface is 
a null vector. Let us consider the Killing vector that generates time translations, ^ = ^^e^. In 
the Schwarzschild spacetime this vector is simply ^ = 64. This Killing vector is normal to the 
horizon, so that S,^^^ = and this is why the event horizon is called a Killing horizon. More 
specifically, ^^S,fi is constant on the horizon thus the gradient V'^(^'^^^) is also normal to the 
horizon. Hence, there exists a function k such that 

^"{e^f.) = -2Ke- (D.l) 

Since the field .^ is a Killing vector it satisfies V^^u + V^^^ = 0, then the above equation can be 
rewritten as 

rv^e. = -i^c^- (D.2) 

By using the hypersurface orthogonal vector property C[^VuS.p] = 0, one can find the relation 

'^' = -^(V^e.)(V^r) (D.3) 

evaluated at the horizon. The Schwarzschild metric is diagonal so we have ^^ = 5^" and ^^ = 
Sfitdtt and for evaluating the covariant derivative, the only nonzero ^u,fM is ^t,r = gtt,r since the 
metric components dependent only on r. Thus we get [187 



-V9«(5«,.)2 = J^=L. (D.4) 

4 2 ^/-gu9rr 



37 



D.2 Conical Singularity 



The expression (D.4) can be obtained from another method. Notice that the horizon property 



imphes that near r = ru, the matric takes the fohowing form 

fiv 

ds^ ~ -A{r - rH)dt^ + t^t r + r^dO^^^ (D.5) 

This metric has a coordinate singularity at the horizon, but only if the Euclidean time t = it 
is periodic with a particular period t r^ t + (3, which one then define the inverse temperature 
(3 = 1/T. The metric read in Euclidean signiture 

dr^ 

ds^ ~ A{r - rH)dT'^ + ^, r + r^dO^ . (D.6) 

B[r - ru) 



With 



A 2 rf^' 2 

"/O =^7 ^' P= —^yf -rn- (D-7) 



the metric (D.6) becomes 

ds^ ~ p^K^dr"^ + dp^ + r^dnl_^, (D.8) 

where k = vAB/2. This geometry is regular if the Euclidean time r is periodic with a period 

/3 = ^. (D.9) 

Here we have defined k and by the first law of the (stationary) black hole thermodynamics, this 
K corresponds to the surface gravity and the black hole temperature is given by 

This will give the same result as the one obtained by direct calculation of the surface gravity 



(D.4), but this Euclidean signature formulation has the advantage of showing subtleties for the 
noncanonical horizon cases, and verifying the relation between Hawking temperature and surface 
gravity. 

E Black hole temperature and entropy 

E.l AdS black holes 

The AdS^ black hole is 

ds^ = ^(-fdt^ + dx^ + ^] (E.l) 

where f{z) = 1 — z^/z^. We expand the denominator of gzz near z ^ z^ as 

Z^f{z) ~ zJrJiZm) + {2Zmf{Zm) + Z^f {Zm)) {z - Zm) 

= Zmf'{Zm){z-Zm), (E.2) 

and for gu part, we have 

/ J [Z-m) J [^mj^rn ~ '^/v^mj^rra / n 

l2 - ~^ ^ li {Z - Zm) 

= ^-^{z-zm). (E.3) 
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Thus the near horizon (Euclidean) metric is 



ds^ ^ L^l^i^{z - Zm)dT' + ^dx'' + 



.2 , ^ j-2 



L' 



7^ 



7^ 



^mJ \^m)\Z Zra) 



dz^ 



If we take 



dp 



L 



dz 



2L 



then, the near horizon metric becomes 



d {^z - Zm) 



l2 

ds ^ K p dr + dp -\ — T^dx 



where 



\nzm)\ 



This also can be calculated by the relation (D.4) directly from (E.l) as 

|2^/-^V'I \f\Zm 



r 1 \dzgtt\ 

K = lim -^^=^= 

z-^zm 2 ^J-guOz 



lim 

Z ^Zm 



2z2 



Now, the temperature is given by 



T 



The area of the horizon is 



A 



I 



27r TXZr. 



73^ 



1^ xJgi. 



L\ 



-Vs 



(E.4) 
(E.5) 

(E.6) 
(E.7) 

(E.8) 

(E.9) 
(E.IO) 



fixed ^m 

where g± is the determinant of the transverse part of the metric and V3 is the volume. Then by 
the first law of the black hole thermodynamics, the entropy is 



S 



A 



L'V3 



4G5 AG^zl' 



(E.ll) 



E.2 Another ansatz 



Now, we consider the following metric ansatz 188 



ds^ = e^^ {-hdt^ + dr) + 



^2^ , ^2B 



dr"^ 
IT' 



fE.12) 



where A, B and h are some functions of r. We assume that the geometry is asymptotically AdS. 
A regular horizon r = th arises when h has a simple zero. It is also assumed that A{r) and 
B{r) are finite and regular functions at r = rn- We consider the metric near the horizon. For 

J- / 5r-r ) 

h{r) _ h{rH) ( h'jrn) _ 2h{rH)B '{rHy 
Since B is finite and hirn) = 0, we obtain 






(E.13) 



h{T) ^ h'JTH) _ 

(r) - p2B(r„)^'^ ^">- 



o2B 



Similarly, for gu part we get 



/i(r)e2^W ~ h'{rH)e^^^''"\r - th). 



(E.14) 



(E.15) 
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Then, the near horizon geometry with Euchdean time t = it becomes 

h'[rH){r-rH) 
With 

dp = ^— = -^jj^=d (V^^^^) , (E.17) 

we arrive at 



Thus, we obtain the surface gravity 

'^ ^''"Ve^^^^^-^^''^^ (E.19) 



Again, using (D.4), this process also can be done directly from (E.12) as 



,. 1 Idrgttl ,. 1 

K = nm ^^^^^ = iim - 

r^TH 2 y/-gttgrr '^^''H 2 

Now we define the temperature 



dr 



Ve2^e2^ 2 



r=^ ^U!:^^ (g_21) 

47r 

Now, the area is j4 = e^M'^n) fQp ^j^j^^ volume and Gat = kI/Stt. Thus the entropy density takes 
the form 

5 = -^ = ?5e3^(^«). (E.22) 

F Hawking-Page transition and deconfinement 

The partition function for canonical ensemble is given by 

^^e-^(0.O^ (F.l) 

where (pel is a classical solution of the equation of motion with suitable boundary conditions. 
When there are multiple classical configurations, we should sum over the all contributions or may 
take the absolute minimum which globally minimizes IsuGRA and dominates the path integral 
for the leading contribution. If there are two or more solutions to minimize IsuGRA-, there 



may be a phase transition between them. This is the Hawking-Page transition 112] , a thermal 
phase transition to a black hole geometry in asymptotically AdS space. Through the AdS/CFT 
correspondence, it was generalized in [3l[l2] that the corresponding dual of this transition is 
confinement-deconfinement transition in boundary gauge theories. 

F.l Hawking-Page phase transition 

The gravity action we study here is 



^ = -16k/^'^^('^+§)- <"-^' 



where G5 is the five-dimensional Newton's constant and R is the radius of AdS^. For a solution 
of the equation of motion, 

Tlpiu - ijTlgf,u - -j^gt^iy = 0, (F.3) 
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the action (F.2) becomes 



' ''d'xV-9=:,-^,V{e) (F.4) 



ds^=[^ + ^] (^^^ + ( 1 + 152 ) ^^^ + ^^^^3' (F-5) 



where V{€) is the volume of the space time. 

There is a solution into AdS space which minimize (F.4). In this case, the AdS^ metric is 
given by 

where dfig is a metric on the three sphere S^. We call this solution as Xi. The Euclidean time 
direction is chosen that the asymptotic boundary at r — )• oo becomes R'^ x S^. Thus the dual 
gauge theory lives on the spatial manifold S''^. 

There is another solution, X2, on a space with a certain temperature. The geometry is now 
the AdS-Schwarzschild black hole and the metric is given by 

ds'={l + i- ^) dr'+{^ + i- ^) " dr^ + r'd^l (F.6) 

where /i = 16ttG5M/3V{S^). The radial direction is restricted to r > r+, where r+ is the largest 
root of the equation 

1 + ^-^=0. (F.7) 

The metric is smooth and complete if r is periodic with the period 

The topology of X2 R^ x S^ and the boundary is R^ x S^. 

The geometry S^ x S^ for at large r from both Xi and X2 configurations are explained. Let 
the 5^ radius be /? = (r/i?)/3o and the S^ radius be /? = r/R. If we wish to make the top ology 



S"^ X S^, we should take /?//? = /3o — )• and this is the limit of large temperatures. From (F.8) 
it seems that this can be done with either r_|_ — )• or r+ — )• 00, but the r^ — )• branch is 
thermodynamically excluded |112j and we have the large r+ branch, corresponding to large M. 
Therefore for either Xi or X2, the topology becomes 5"^ x S^ at large r. 

Actually, both I{Xi) and I{X2) are infinite, so we compute 1(^2) — -^(^1) to get a finite 
result. We put a cut-off Rq in the radial direction r. Then, the regularized volume of the AdS 
spacetime for Xi is given by 

rl3' rRo r 

Vi{Ro) = dr dr dQr^. (F.9) 

Jo Jo Js^ 

And for X2 case of the AdS-Schwarzschild black hole, it is 

/■fl) fRo f 
V2{Ro) = dr dr dQr^. (F.IO) 

Jo Jr+ Js3 

In order to compare the I{Xi) and /(X2), we match period of r so that the proper circumference 
of the Euclidean time direction at r = Rq is the same each other. This can be done by setting 



/3'V^ + l = /3oV^ + l-^ (F.ll) 

at r = i?o and we determine (3' . Then the difference I{X2) — liXi) becomes 

lim {V2 - Vi) = ^^ +,^ — ^. (F.12) 



2ttG5R^ Ro^oo^ ^ 4G5(4r2 +2i?2) 



Then (F.12) changes its sign at r-|_ = i? and the phase transition to the AdS-Schwarzschild black 



hole geometry takes place. 
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r=r, 




(a) thermal AdS (b) AdS-Schwarzschild 

Figure 7: The topology of the thermal AdS geometry and the AdS-Schwarzschild geometry. 

F.2 Confinement and deconfinement 

Deconfinement at high temperature can be understood by the spontaneous breaking of the center 
of the gauge group. The corresponding order parameter is the Polyakov loop that is defined by 
a Wilson loop wrapping around r direction as 



Nc 



(F.13) 



where P denotes the path ordered configuration. In our case the gauge group is SU{Nc) and its 
center is I^Nc- 



The gravity dual calculation for the expectation value (P) was performed in [189 190 from 
the regularized area of the minimal surface ending on the loop, 



{V) 



-fiA 



(F.14) 



where ^ is the fundamental string tension, and A is the area of the minimal surface ending on 
the loop. The boundary of the spaces Xi and X2 is S^ x S^ and r direction wraps around 
the loop C. To compute (V^C)), we should evaluate the partition function of strings with its 
worldsheet D which is bounded by the loop C. 

For the low-temperature phase, the space is Xi and the geometry is the thermal AdS with 
topology S^ X B^. Then the loop C is never contracted to zero in Xi and C is not a boundary 

This immediately implies that 



of any string worldsheet D, see Figure 7(a 



{r) = o 



(F.15) 



in the large Nc limit. Thus at the low-temperature, the thermal AdS geometry is stable and the 
configuration Xi corresponds to the confined phase. 

On the other hand, for the high-temperature phase, the relevant space is X2 and the geometry 
is the AdS-Schwarzschild black hole with topology R^ x S^. In this case, the loop C can be a 



boundary of a string worldsheet D = B^, see Figure 7(b) Then the (regularized) area of the 
surface gives the non-zero Polyakov loop expectation value. 



(V) / 0. 



(F.16) 



Therefore, at the high-temperature, the the AdS-Schwarzschild black hole geometry is more 
stable and the configuration X2 corresponds to the deconfined phase. 
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